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DECUS Program Library Wrife-up 


DECUS NO. 8-652 


REGRESSION ANALYSIS PACKAGE 

FOR THE POP 8 WITH 4K OF CORE AND AN A5R33 CONSOLE 

In Numerical Analysis , we have a therom stating that 
any curve may be represented by a polynomial equation ; if, 
you have enough terms in the equation. In actual practice, we 
can fit most curves with only a few terms , In this program, 
we assume that a dependent variable W may be represented by 
a polynomial function of independent variables ( X, Y, Z), 
We enter data for many points , and ask the computer to find 
the coefficients that will give a least squares fit , We 
specify the number of terms in the equation , This we call 
THE DEGREES OF FREEDOM. 

This package is a group of programs for making a multiple 
regression with up to 3 independent variables having up to 
20 degrees of freedom. They are useful in two ways: 

1, If you want the computer to find a formula that will 
calculate data given by a family of curves . For example, 
you might want to run a heat balance on a boiler . You 
would want the computer to determine the enthalpy of steam 
for many different values of temperature and pressure . Also, 
you might want to read the guarantee curves supplied by the 
boiler manufacturer. 

2 . You might want to find out if there is a pattern in a 
large mass of test data . These regression programs will 
find the pattern ( if it exists ) and will print the 
standard deviation from the pattern. 

These regression programs are delivered on four tapes. 

RP The first is the RESIDENT PACKAGE that we keep in core 

5000- at all times for use also on other programs . It includes 

6777 the 3 Page — 27 Bit Floating Point Package DECUS B-375B. 

RA1 The second, -- REGRESSION ANALYSIS part 1 — calculates the 

7000- regression coefficients. 

7577 

RA2 The third -- REGRESSION ANALYSIS part 2 -- will reread 

4400- the data and expand the polynomial equations using the 

4777 regression coefficients of part 1 . The calculated value 

of the dependent variable V is compared with the given 
value of W o 

QD&flL The fourth QUICK DUMP and QUICK LOAD- are in RIM format 

7600- as they overwrite the BIN loader . Since the regression 

7755 parameters must be loaded in a relocatable format , the BIN 

loader could not be used. 

On the next page, we explain the regression formulas . We 
follow with a list of the programs in the order that they 
appear in core. We will discuss them in this order. 











REGHECSIOrJ 


EQUATIONS 
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WE WILL SIMPLIFY THIS EXPLAI'IATIOM BY USING A SPECIFIC 
FOP^MULA RATHER THAN A GENERAL ONE. WE WILL USE THE UP ARROW 
(t) TO DENOTE EXPONENTIATION, AND THE PRIME <*) TO DENOTE 
DIFFERH'JTIATION. V IS A POLYNOMIAL FUNCTION OF (X & Y) . VARIABLES 
A, B, C, D, E, & F ARE THE CO EFFI Cl HMTS TO BE MULTIPLIED BY EACH 
SUCCESSIVE TERM. 

V = A + BX +CXt2 +DY +EYX +FYXt2 

IN THE ABOVE EQUATION, V HAS 3 DEGREES OF FREEDOM 
WITH RESPECT TO X < A PARABOLA ) AI'JD 2 DEGREES OF FREEDOM 
WITH RESPECT TO Y (A STRAIGHT LINE ) . WE WANT TO FIND VALUES 
FOR A, E, C, D, E, £ F THAT WILL GIVE A MINIMUM ERROR (SQUARED) 
UHe\I SOLVED WITH AI'JY GIVEN VA1.UES FOR X AMD Y. W WILL BE 
A GIVEN OR MEASURED VALUE OF THE DEPENDENT VARIABLE. V IS 
THE CALCULAiTED VALUE GIVEN BY THE ABOVE FORMULA. 

THE ERROR SQUARED IS: U = CV-W)t2 

FOR A MINIMUM VALUE OF U , THE DERIVATIVE WITH RESPECT TO A 
WILL BE ZERO. 

UVA’ = 22(V-W) WA', BUT VVA* = 1 ,* SO XV = XW 
AND 'SA + 5BX +XCXt2 +'2DY +2EYX +rFYXT2 = 2W 

IN A LIKE MANNER, WE CAN DIFFERENTIATE WITH RESPECT TO B AND 
GET ANOTHER EQUATION : 

UVB' = 2S.(V-V) VVB* ; BUT VVB* = X J SO , 2VX = :2WX 

XAX +'2EXt 2 +XCXt 3 +SDYX +2:EYXt 2 +SFYXt 3 = 2WX 

WE CAN CONTINUE THIS PROCESS AND GET 6 EQUATIONS WITH 6 
LNKNOWI\I COEFFICIENTS. WE WILL ARRAiNGE THESE EQUATIONS IN 
MATRIX FORM. 


REGRESSION MATRIX 


N 

+2X 

+•50 (t 2 

+5.Y 

+5YX 

+5YXt 2 

=s:w 

?x 

+FXt 2 

+XXt 3 

+SYX 

+5YXt 2 

+5YXt 3 

It 

2Xt 2 

+SXt 3 

+5Xt A 

+5YXt 2 

+5.YXt 3 

+5YXt A 

=5:wxt 2 

5Y 

+XYX 

+XYXT2 

+¥Yt2 

+5Yt 2X 

+r'i't2XT2 

=SV7Y 

5YX 

+XYXt 2 

+-SYXt 3 

+-SYt 2X 

+5Yt 2Xt 2 

+5Yt2Xt3 

= Z\'PCX 

lYXt 2 

+S:YtXt 3 

+5YXt A 

+5Yt2Xt2 

+'5Yt 2XT 3 

+2Yt 2XtA 

=-S:WTXt 2 
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LI5T OF PR0GRAM5 

4400—4450 

5-4420 

DRC 

DUMP REGRE55I0N C0EFFICIENT5 to tape. 

4451 —4 577 
5-4472 

R;A2 

REGRE55I0N ANALY5I5 part 2. This program rereads 
the data and calculates the error for each data 
point. It prints the RM5 error or the 5TANDARD 
DEVIATION at the end. 

4600—4777 
JM5 I A4600 
PLY 

POLYNOMIAL EXPAN5I0N.' This routine will calculate 
the dependent variable, by the polynomial expansion 
of the regression coefficients. 

5000—5021 
JM5 I A5000 
ME5 

ME55AGE. This routine may be used to type a 
message that was loaded in core by QL at the 
address in the AC. 

5041—51 77 
CR, POA, 

FPT 

CARRIAGE RETURN, PRINT OCTAL ACCUMULATOR, AND 
FLOATING POINT TRACE. These are utility programs. 

5200—5227 

5-5200 

PME 

POOR MAN'5 EDITOR to copy a tape with some control 
over the typewriter. 

5233 — 5337 
JM5 I A5271 
ANIL 

ANTILOG. This program will calculate the anti¬ 
logarithm of the Naperial logarithm given in the 
floating point accumulator. 

5400—6577 
JM5 I A5400 
FPNT 

Five page FLOATING POINT PACKAGE with floating 
output — DECU5 8-375B. 

6502—6565 
JM5 I A6502 
RINT 

READ INTEGER. This program reads integer numbers, 
from the keyboard. 

6600—6660 
JM5 I A6600 
FIN 

FLOATING INPUT routine to read floating point 
numbers from the keyboard. 

6661 —6777 
JM5 I A6661 
LOG 

Routine to calculate the Naperian LOGARITHM of 
the number in the floating point accumulator. 

7000—71 77 
JM5 I A7012 
55E 

50LVE 5IMULTANE0U5 EQUATIONS. This routine will 
solve up to 28 simultaneous equations. 

7200—7577 

5-7342 

RA1 

REGRESSION ANALYSIS part 1. This program builds 
the regression matrix, and calls 55E to solve 
for the regression coefficients. 
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7600—7675 

5-7600 

QD 


7700—7753 

5-7700 

QL 


LI5T OF PR0GRAM5 continued 


QUICK DUMP. This program will dump one page of 
core starting at the address in the switch register. 
It will continue on the next page every time you 
depress CONT after each halt. Two characters will 
be typed to define each word. Each line will start 
with a bell (207) followed by the address at which 
the line starts in core. 

QUICK LOAD. This program will load a tape cut by 
QD or by DRC. You must type the loading address 
before loading a tape cut by DRC. ( if you forget, 
you may clobber core. ) 

The two proceeding programs are on a tape in RIM 
format for loading by the R:IM loader at 7756, 
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4400-4450 DRC - DUMP REGRESSION COEFFICIENTS —S-4400 

This program is automatically entered on exit from RA2. 

A halt at the beginning of the program will give you time to 
turn the punch on. If you know where the coefficients are to 
be stored, you should type BELL (207) ADDR before you depress 
CONT. This will tell the loader QL where to store the coefficents. 
You should also set some identifying pattern in the switch register. 
This pattern will be stored in N4 where it may later be examined 
to identify the data. All other parameters are set by the program 
from data entered in RA1.The following parameters are put on tape, 

NL * a parameter to tell PLY whether or not to take the 

logarithm of a variable before calculating the dependent 
variable. 




NL = 0 1 

2 

3 

4 


Take 

antilog of W ? 

YES. YES 

YES 

YES NO 


Take 

log of 

X ? 

YES YES 

YES 

NO 

NO 


Take 

log of 

Y ? 

YES YES 

NO 

NO 

NO 


Take 

log of 

z t 

YES NO 

NO 

NO 

NO 


N1 = 

Negative 

number 

of degrees of 

f reedom 

for variable 

X . 

N2 = 

Negative 

number 

of degrees of 

freedom 

for variable 

Y . 

N3 = 

Negative 

number 

of degrees of 

freedom 

for variable 

Z . 

N4 = 

Identification. 

Arbitrary pattern 

set 

in switch register. 

N5 = 

Number of 

words 

required to store 

the 

regression coefficients. 


The regression coefficients follow. 


4451 -4577 -RA2 — REGRESSION ANALYSIS part 2 — S-4472 

This program will reread all of the input data. For 
each point, a value V will be calculated using the regression 
equation with the coefficients determined in part 1 .See page 2 . 
The error (V-W) will be determined. For each point, the following 
values will be printed as shown; 

V (V-W) POINT NO. 

W X Y Z 


At the end of the run, the RMS error will be typed. This program 
may be entered by depressing CONT after the halt at the end 
of part 1 at 7264. However, if you specified more than 26 
degrees of freedom, this program will have been overlayed by 
the matrix built in part 1. So, you must reload this program 
and start at 4472. 







4600-4777 PLY POLYNOMIAL EXPANSION 


JM5 I A4600 
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This routine will expand the polynomial equation. On return 
the floating point accumulator will hold the calculated value V 
You will need this routine in your final program that is to read 
the data correlated by the regression coefficients. Use the 
following sequence of instructions in your program. 


ANL, 

NL-1 


/ ADDRESS OF NL-1 

A4600, 

4600 


/ ADDRESS OF PLY 

AOX, 

X 


/ ADDRESS OF X 


CLA 




TAD 

ANL 

/ GET ADDRESS OF NL-1 


JMS I 

A4600 



AOX 


/ ADDRESS OF X 




/ CONTINUE WITH PROGRAM. 




/ FAC WILL CONTAIN V 

X, 

0000 


/ INDEPENDENT VARIABLES 


0000 


/ LOCATED FROM THIS POINT ON. 


0000 



Y, 

0000 




0000 




0000 



z. 

0000 




0000 




0000 



NL, 

0000 


/ REGRESSION PARAMETERS FOLLOWED 

N1 , 

0000 


/ BY REGRESSION COEFFICIENTS 

N2, 

0000 


/ LOCATED FROM HERE ON IN ORDER. 

N3, 

0000 


/ QL WOULD BE USED TO LOAD THE 




/ DRC TAPE TO PLACE THEM HERE. 

5000-5021 

ME 5 

MESSAGE 

— JMS I A5000 


This routine may be used to print a message that was 
loaded in core by QL at the address given in the AC. If 
your message includes : 

1. Back arrow up arrow ( ^ f ) it will generate a 
carriage return line feed when printing, 

2, Two back arrows ( ) will terminate the message 

Each pair of symbols, given above, must be contained in one 
word. When using QUICK LOAD to load your message, the AC 
will contain 7777 if you are entering the first character 
of the pair. Otherwise, it will be blank. So, if the AC 
is blank when you want to enter the first back arrow ( ^ , 

hit a space first. 
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5042-5105 CR CARRIAGE RETURN — JM5 I 74 

This routine will issue a carriage return - line feed. 
The accumulator will be cleared. It will also keep a line 
count and type 6 blank lines whenever 60 lines have been 
counted. If you call this routine with 7777 in the AC, 
it will skip to the next page. 


5106-5133 POA PRINT OCTAL ACCUMULATOR — JMS I 73 


This routine will cause 
to be printed in octal. 


the contents 


of the accumulator 


6305- 

■6312 

TYPE 

— 

JMS I 75 

This 

routine 

will 

cause 

the character in the accumulator to 

be typed and 

the 

accumulator cleared. This routine furnished 

with 

DECUS 

0-375B was changed as follows to improve timing. 

6305 

0000 

TYPE, 

0 


6306 

6041 


TSF 


6307 

5306 


JMP 

.-1 

6310 

6046 


TLS 


631 1 

7200 


CLA 


6312 

5705 


JMP I 

TYPE 

5134-5177 

FPT 

FLOATING POINT TRACE — S-5134 


To use this trace, you must first establish the linkage. 
You can use ODT ( Octal Debugging Technique ) to insert codes 
as follows. 


1 . 

Insert 

51 34 

at 

address 

0077 

2. 

insert 

5477 

at 

address 

5404 

3. 

insert 

LOW 

at 

address 

5050 

4. 

Insert 

HIGH 

at 

address 

5051 


LOW is the lowest address to be traced. 
HIGH is the higest address to be traced. 


All floating point interpretive instructions that fall in 
the range between LOW and HIGH will be traced. First, the 
carriage will be returned, then the address will be typed in 
octal. Finally, the FAC will be typed in E format. 

LOW and HIGH will usually be on the same page. But, 
this is not necessary if they are both on the same side of 
address 4000. We cannot have one positive and one negative 
number used to specify the tracing range. On the next page, 
we show you two traces that were run on the ANTL routine 
and one on the LOG routine. 






■5 


B 


(-. 068 ) 

1620G trace e 
-.868 A 


5273 - 1.252259E+00 
521ii - 1.252259E+00 

5303 2.52259 1E-01 

5304 2.52259 1E-01 

5306 2.52259 1E-ei 

5307 1.000000E+00 

5310 7.477409E-01 

531 1 7.477409E-01 

5324 2.500000E-01 

5325 7.477409E-01 

5331 2.591472E-01 

5332 2. 59 1472E-01 

5333 6.715728E-02 

5334 6.015735E+01 

5335 6.015735E+01 

5336 -6. 018043E+02 

5337 - 1.000384E+01 

5340 - 1.026298E+01 

5341 1.752033E+00 

5342 1.752033E+00 

5343 2.59 1472E-01 

5344 5. 182944E-01 

5345 2.958246E-01 

5346 1.295825E+00 

5347 1.295825E+00 

5350 1.679162E+00 

5351 4. 197904E-01 


1620G 
64 L 
6663 
6707 
6710 

6716 

6717 

6720 

6721 

6722 

6723 

6724 

6725 

6726 

6727 

6730 

6731 

6732 

6733 

6734 

6735 

6736 

6737 


4. 197904E-01 


(2.00148) 

2.00148 A trace e' 

5273 2.887525E+00 

5274 2.887525E+00 

5303 -8.875249E-01 

5304 -8.875249E-01 

5305 -8.875249E-01 

5324 4.000000E+00 

5325 -8.875249E-01 

5331 3.075927E-01 

5332 3* 075927E-01 

5333 9.461326E-02 

5334 6.018480E+01 

5335 6.018480E+01 

5336 -6. 018043E+02 

5337 -9.999273E+00 

5340 - 1.030687E+01 

5341 1.708151E+00 

5342 1.708151E+00 

5343 3.075927E-01 

5344 6.151854E-01 

5345 3.601469E-01 

5346 1.360147E+00 

5347 1.360147E+00 

5350 1.850000E+00 

5351 7. 399999E+00 7.399999E+00 


trace log(64) 

6.400000E+01 
7. 000000E+00 
6.400000E+01 
5. 000000E-01 
1.207107E+00 
1.207107E+00 
5.000000E-01 
2.071068E-01 
1.715729E-01 
1.715729E-01 
2V943725E-02 
2.943725E-02 
5.989786E-01 
1. 763228E-02 
9.79 1028E-01 
2.882209E-02 
2.914213E+00 
4.999999E-01 
9.999999E-01 
6.000000E+00 

4. 158883E+00 4. 158883E+00 









5200-5227 


POOR MAN'S EDITOR 


5-5200 
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PME 


This is a short program that will copy a tape with some 

control over the typewriter. 

1, if you set a character code in the switch register, the 
typewriter will halt before typing the selected character. 
If you want to type the selected character, depress CONT. 
If you want not to type the character, depress EXAM, CONT. 

2, If you put SWO up, every character will be selected; so, 
you can step accross the page — one character at a time, 
each time you depress CONT, 

3, If you want to add a few characters to the output tape 
that were not on the input tape, you can turn the reader 
off and depress CONT. All characters typed will appear 
on the output tape. 

4, If you put SW3 up, it will temporarily cancel your 
character selection, and the typewriter will keep going. 
When you get near the point where you want to stop; put 
it down and the typewriter will halt before typing the 
selected character. 


5234-5270 FIX — JM5 I A5234 

This routine will shift the FAC left or right untill the 
decimal point registers at the low order end of AC1. AC2 
and AC3 will be set to zero. Any digits shifted out of the 
high order side of AC1 will be lost and no error flagged. The 
routine will return with the integer portion of the number in 
the machine accumulator. If the sign was negative, it will be 
two's complimented. 


5271-5377 ANTE ANTILOGARITHM — JMS I A5271 

On return from this routine, the floating point accumulator 
will be replaced with the antilogarithm {e ) of the Naperian 
logarithm that was in the FAC when called. 

6502-6565 R'lNT READ INTEGER — JMS I A6502 

This routine is used by FIN below to read in the decimal 
digits and convert to binary. Therefore, it uses the conventions 
described below. It will load the number in AC3 overflowing in¬ 
to AC2, and thence into AC1. If you keep your integers smaller 
than 2040 , you can use single precision. You should test SGN 
at addr. 51 and compliment if it is negative. 




6600-6660 


FIN 


FLOATING INPUT 


JM5 I A6600 


This routine will read decimal digits from the keyboard, 
convert to binary, and load the number in the FAC. E FORMAT 
IS NOT AN ACCEPTABLE INPUT FORM. 

Too many digits at input will produce an unflagged error. 

We cannot handle 11 digits. 

34567890120. I 2.081518E+08 Unflagged error 

unless the first digit is 2 or less 

23456789010. I 2.345679E+10 

29876543210 I 2.987654E+10 


Leading zeros don't count 
00000987654321 I 9.876543E+08 

.9876543212345 I 1.529630E-03 Unflagged error 
.0000098765432 I 9.876542E-06 



The sign may appear at either end of the number. 
-123.456 is the same as 123,456- 


Other than minus sign (-) and decimal point (.) , all 

characters coded less than (260) are ignored. Carriage return 
(215), line feed (212), plus sign (253), and space (240) are 
all ignored by the input routine. 


Any character coded greater than nine (271>) will terminate 
the input. The terminating character will be saved in address (47) 
Your program should examine this register to identify whether 
and where the number is to be stored. Thus, you can arrange to 
have incorrect numbers overwritten by corrected ones. 



6661-6677 LOG LOGARITHM —JM5 I A6661 


On return from this routine, the floating paint accumulator 
FAC will be replaced with the Naperian logarithm of the number 
that was in the accumulator when called. A halt at 6667 will 
occur if the FAC is zero or negative when called. 

The logarithm and antilog routines give about siz digit accuracy. 


LOGARITHMS 


ANTTLOGARITHMS 


.00000987654321 

23456789010. 

.4197904 

7.4 

987654321 

98765432 


L- 1. 152535E+01 
L 2.387843E+01 
L-8.679998E-01 
L 2.001480E+00 
L 2. 071084E+01 
L 1.840826E+01 


-11.52535 
23.87843 
-.868 
2. 00148 
20.71084 
18.40826 


A 9.876536E-06 
A 2.345683E+10 
A 4. 197904E-01 
A 7.399999E+00 
A 9.876489E+08 
A 9.87654 IE+07 



7000-7177 SSE 


SOLVE SIMULTANEOUS EQUATIONS — JMS I A7102 


1 1 


This routine will solve a set of simultaneous equations 
that have previously been loaded in matrix form starring at 
address 236. The accumulator must be clear at entry. The 
single precision word at address 100 must be set equal to 
the negative number of equations. The word at address 105 must 
be a positive number equal to the number of equations,multiplied 
by 3 , It shows how many words are needed to store one column 
of the matrix. The unknowns ( or answers) will be printed before 
exit. Also, the answers will be stored as floating point numbers 
( tripple precision ) starting at address 106. The matrix of 
equations must be stored in column major order ( as in Fortran) 
with the given constants to the right of the equal sign stored 
in the last column. 


7200-7377 RAI REGRESSION ANALYSIS part 1 S-7342 

This program will first read four parameters —A, B, C, &. D, 

A = number of degrees of freedom for variable X 

B = number of degrees of freedom for variable Y 

C = number of degrees of freedom for variable Z 

D will tell which variables are to be converted to logarithms 
before use. 





D 


0 

1 

2 

3 

4 

T eke 

log 

of 

W 

7 

NO 

YES 

YES 

YES 

YES 

Take 

log 

of 

X 

? 

NO 

NO 

YES 

YES 

YES 

Take 

log 

of 

Y 

7 

NO 

NO 

NO 

YES 

YES 

T ak e 

log 

of 

Z 

? 

NO 

NO 

NO 

NO 

YES 


Some data will correlate better if plotted 
or semilogarithmic graph paper. For this kind of 
set D at a value that will cause the logarithm 


on logarithmic 
data, you should 
to be taken. 


These parameters may be read in any order. If you find 
that you have entered a wrong number, you may overwrite it 
with a corrected one. When you have all parameters correctly 
entered, enter an E to terminate the list. 


To enter the above parameters, you simply type a number 
followed by the designating letter A, B, C, or D. Space, Carriage 
Return and Line Feed will be ignored. If not entered. A, B, &. C 
will default at (1.) . D will default at zero. 

Now that the parameters are entered, you may enter the 
dependent variable W; and the independent variables X, Y, &. Z. 

For each data point, a value for W and one for X should be 
entered; with Y and Z optional. 
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Here also, you may overwrite any variable incorrectly 
entered with a corrected number. Variables may be entered in any 
order. You type the number, followed by the designating 
character Wy X, Y, or Z, Space, Carriage Return , of Line Feed 
will be ignored. The data for one point is terminated with 
a left bracket (IT ). If data for a variable was not entered, 
the previous value will be u&ed. Any character coded less than 
(327) will cause its associated number to be discarded. You 
may terminate the run with a back slash (\ ) • 

The program assumes that the dependent variable Vf may 
be repxesHTvted as a polynomial function of the independent 
variables X, Y, &. Z .It will develop regression coefficients 
for use in the equation as explained on page 2 , It will 
develop coefficients to give a least squares error. 

You must enter more data points than there are degrees 
of freedom. In other words, the number of data points must be 
more than the product (A x B x C ). You ought to have many 
more. If you do not have enough points, you may halt on a 
divide by zero at addr. 6064, 

The example in the appendix will help to explain the 
entry of data. You will probably want to cut a tape in the 
local mode for use in both part one and part two. You can 
use the POOR MAN'S EDITOR — PME at 5200 to make minor 
corrections to your data tape. 

This program calls REG to form the matrix, and* S5E 
to solve for the regression coefficients. 


7413-7577 REG REGRESSION MATRIX — JMS I A7413 

This routine will read the values for the dependent 
variable W and the independent variables ( X, Y, &. Z ) 
and form the regression matrix as shown on page 2 . It will 

return to RAI after the endfile character — back slash { \ )—— 
is read , 



The next two programs are on a tape in RIM format. 
They will overlay the BIN loader ; so they must be 
loaded by the boot at 7756. They are not really part of 
the REGRESSION ANALYSIS SYSTEM , but are needed to put 
the regression coefficients in core at relocatable 
addresses for use by other programs. 


7600-7677 QD QUICK DUMP — 5-7600 

You can use this program to dump or display any part 
of core memory at 2 characters per word. This is useful 
when you are debugging a program that changes its own 
instructions, and you want to find out where it went wrong. 
You take a listing of any page in core in about 30 seconds. 
You can compare it with a previous listing before it went 
wrong. 


To start the dump: 

1. Put 7600 in memory address register. 

2. Set the first address to be displayed in the switch reg. 

3. Depress CONT. 

You can use the table on the next page to translate any 
pair of characters in the listing into the equivalent octal 
word. 

The program will halt after 120 words have been listed. 
Depress CONT to go on to the next page. Each line in the 
listing will be preceeded by BELL and a two character address 
where the line is to be loaded. So, if a tape is cut during 
the dump, it may later be reloaded by QUICK LOAD. 

7700-7753 QL QUICK LOAD — 5-7700 

This program may be used to reload any part of memory 
that was dumped by QUICK DUMP or DRC. The tape dumped 
by DRC will not contain a starting address; so, vou 
must supply it. Before starting'the reader, type BELL followed 
by the 2 chiaracter address where'you want to start - loading. 
Using the translate table on the next page,,you can translate 
the octal address into its 2 character equivalent. 

Each character loaded will be echoed on the typewriter. 
But, if you want to mute the echo, you type: 

BELL 

before starting the reader. After the tape has been loaded, 
you can restore the echo by typing; 


BELL >^8B: 





u 


SYMBOL TABLE 


IN THIS TABLE, EACH OCTAL INSTRUCTION IS FOLLOWED BY ITS EQUIVALENT 
CHARACTER PAIR, AND THEN BY ITS SIX CHARACTER NUMONIC EQUIVALENT. 

NOTEi & - AIv’D ; V = OR 


7100 

Y 

CLL 

7240 

Z« 

CLACMA 

7450 

\H 

SNA 

7550 

3H 

SPA&NA 

7000 

X 

NOP| 

7101 

Y! 

CLLI AC 

7300 

C 

CLL CL A 

7460 

\P 

SZAVNL 

7600 

t 

CLA 

7001 

X! 

lACj 

7104 

Y$ 

CLLRAL 

7401 

\l 

NOP 

7470 

\X 

SNA&ZL 

7604 

t$ 

LAS 

7002 

X” 

b:» 

7106 

Y& 

CLLRTL 

7402 

\” 

HLT 

7500 

] 

SMA 

7621 

t 1 

CAM 

7004 

x$ 


7110 

Y< 

CLLRAR 

7404 

\S 

OSR 

7501 

1 ! 

MQA 

7640 

t0 

SZACLA 

7006 

x& 

RTl. 

7112 

Y* 

CLLRTR 

7410 

\< 

SKP 

7510 

] < 

SPA 

7650 

tH 

SNACLA 

7010 

x< 

RAR 

7120 

Y0 

STL 

7420 

\0 

SNL 

7520 

]0 

SMAVNL 

7700 


SMACLA 

7012 

X=*: 

RTR 

7200 

Z 

CLA 

7421 

\1 

MQL 

7521 

1 1 

SWP 

7701 

- ! 

CLAMQA 

7020 

X0 

CML 

7201 

Z! 

CLAIAC 

7430 

\8 

SZL 

7530 

38 

SPA&ZL 

7710 


SPACLA 

7040 

x§ 

CMA 

7204 

ZS 

GLK 

7440 

\§ 

SZA 

7540 

] 0 

SMAVZA 

7721 

-1 

CLASWP 

7041 

XA 

CIA 


IN THE FOLLOWING TRANSLATE TABLE, AN OCTAL PAIR IS PLACED DIRECTLY OVER 

ITS EQUIVALENT CHARACTER 


& AND 
^IRIND 

f TAD 
DIRIND 

I ISZ 
DIRIND 

D DCA 
DIRIND 

C JMS 
DIRIND 

J JMP 
DIRIND 



0 1 

01234567 

1 1 

01234567 

<)♦+,-./ 

2 1 

01234567 

01234567 

3 ' 

’01234567 

89:;<=>? 

4 ' 

01234567 

eABCDEFG 

5 ' 

01234567 

IHIJKLMNO 

6 

01234567 

PQRSTUVW 

7 

01234567 

XYZC\]t- 



DECIMAL EQUIVALENT 
OF AN OCTAL DIGIT 


512 

64 

8 

1 

1024 

128 

16 

2 

1536 

192 

24 

3 

2048 

256 

32 

4 

2560 

320 

40 

5 

3072 

384 

48 

6 

3584 

448 

56 

7 
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If you want the tape to halt the computer after loading, 
punch the following characters at the end of the tape: 

BELL-^+X " 


In the above, we used QUICK LOAD to patch itself. 
Of course, any program may be patched in the same way. 
Suppose you want to insert the following patch in your 
program; 

QUICK LOAD 
CHARACTERS 


SYMBOLIC 

LISTING 

DCA N2 
TAD I 11 
JMS I PW 
ISZ N2 


ODT 

LISTING 

5451 /3032 

5452 /I 411 

5453 /4661 

5454 /2032 


BELL 

LI 

8: 

FQ 

0: 


The BELL — control G — introduces the address into 
which the next instruction is to be loaded. The address and 
the instructions that follow must be translated by you from 
octal words to character pairs using the translate table on 
the preceeding page. 
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0200 

/6032 

tsm. 

kcc 

0201 

/6040 


tfl 

0202 

/1206 


tad 

0203 

/4607 


jms 

0204 

/7402 


hit 

0205 

/5202 


jmp 

0206 

/021 0 

ma 9 


0207 

/5000 

ames, 


0210 

/lilt 

cr^ 


0211 

/5657 

no 


0212 

/6700 

M/ 


0213 

/5163 

is 


0214 

/0064 

t 


0215 

/5045 

he 


0216 

/0064 

t 


0217 

/5155 

im 


0220 

/4500 

e 


0221 

/4657 

f o 


0222 

/6200 

r 


0223 

/4154 

al 


0224 

/5400 

1 


0225 

/4757 

go 


0226 

/5744 

Q'd 


0227 

/0055 

m 


0230 

/4556 

en 


0231 

/lilt 

cr* 


0232 

/6457 

to 


0233 

/0043 

c 


0234 

/5755 

om ’ 


0235 

/4500 

B 


0236 

/6457 

to 


0237 

/0064 

t 


0240 

/5045 

he 


0241 

/0041 

a 


0242 

/5144 

id 


0243 

/0057 

o 


0244 

/4600 

f 


0245 

/6450 

th 


0246 

/4551 

ei 


0247 

/6200 

r 


0250 

/6041 

pa 


0251 

/6264 

rt 


0252 

/7100 

y 


0253 

/nil 

eof • 



/ test message 

Flf9 

i ames 

.-3 






DUMP REGRESSION COEFFICIENTS 
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4400 

/0000 

pw, 

0 


4401 

/3047 


dca 

ch 

4402 

/1047 


tad 

ch 

4403 

/7002 


bsw 


4404 

/4210 


jms 

tc 

4405 

/1047 


tad 

ch 

4406 

/4210 


jms 

tc 

4407 

/5600 


jmp 

i pw 

4410 

/0000 

tc , 

0 


Ml 1 

/0250 


and 

p77 

4412 

/ 1215 


tad 

p240 

4413 

/4475 


type 


4414 

/5610 


jmp 

i tc 

4415 

/0240 

p240. 



4416 

/7740 

m40. 



4417 

/7402 


hilt 


4420 

/6032 

drc, 

kcc 


4421 

/6040 


tfl 


4422 

/4474 


cr 


4423 

/ 1216 


tad 

m40 

4424 

/3043 


dca 

int 

4425 

/7604 


las 


4426 

/3104 


dca 

1 04 

4427 

/7327 


7327 


4430 

/I 105 


tad 

1 05 

4431 

/7041 


cia 


4432 

/3046 


dca 

cnt 

4433 

/ 1250 


tad 

p77 

4434 

/3010 


dca 

1 0 

4435 

/ 1410 


tad 

i 10 

4436 

/4200 


jms 

pw 

4437 

/2043 


isz 

int 

4440 

/5244 


jmp 

.+4 

4441 

/4474 


cr 


4442 

/ 1216 


tad 

ra40 

4443 

/3043 


dca 

int 

4444 

/2046 


isr 

cnt 

4445 

/5235 


jmp 

.-10 

4446 

/7402 


hit 


4447 

/5246 

pTT. 

jmp 

.-1 

4450 

/0077 




/ print word 


/ type character 


/ dump regression coefficients. 


/count words / line 
/ set at 32 
/ code from switch reg. 

/ set in 104 for identification 

/ 6 

/' number of coefficients 

/ set count 
/ starting address 


/ 32 words per line 

/carriage return 
/ reset words per line 


/ repeat 

/ halt when done 





MULTIPLE REGRESSION 


Part 2 


4451 

/4474 

exit, 

cr 


/ 

carriage return 

4452 

/4476 


f pnt 




4453 

/0356 


f get 

tct 

/ 

get total error sq-. 

4454 

/4267 


f div 

num 

/ 

calculate average 

4455 

/0000 


f ext 




4456 

/4665 


jms i 

alog 

/ 

calculate square root 

4457 

/7240 


cla cma 



4460 

/ 1061 


tad 

exp 



4461 

/3061 


dca 

exp 



4462 

/4666 


jms i 

aantl 



4463 

/4752 


jms i 

afout 

/ 

print root mean square 

4464 

/5217 


jnip 

dxc-1 

/ 

end of run 

4465 

/6661 

adog, 





4466 

/5271 

aantl 





4467 

/2044 ■ 

num, 





4470 

/0000 






4471 

/0000 






4472 

/7307 

prt2, 

7307 


/ 

-4i part two 

4473 

/ 1 104 


tad 

n4 



4474 

/3100 


dca 

nO 

/ 

set log test parameter 

4475 

/3104 


dca 

n4 

/ 

zero n4 

4476 

/6032 


kcc 




4477 

/6040 


tfl 




4500 

/4476 


f pnt 




4501 

/0754 


f get 

i azro 

/ 

initialize 

4502 

/7356 


f put 

tot 

/ 

total 

4503 

/7773 


f put 

i anum 

/ 

and number 

4504 

/0000 


f ext 




4505 

/4747 


jms i 

aw/xyz 

/ 

read data w, x, y, &. z 

4506 

/5251 


jmp 

BX it 



4507 

/ 1353 


tad 

P77 

/ 

set index register to 

4510 

/4750 


jms i 

aply 

/ 

read pars, in ply 

451 1 

/7567 


adx 


/ 

address of x used by 

4512 

/4474 


cr 


/ 

issue carriage return 

4513 

/4476 


f pnt 




4514 

/7771 


f put 

i awe 

/ 

save calculated value 

4515 

/0755 


f get 

i aone 



4516 

/ 1773 


f add 

i anum 



4517 

/7773 


f put 

i anum 

/ 

increment number 

4520 

/0771 


f get 

i awe 



4521 

/2774 


f sub 

i adw 



4522 

/7772 


f put 

i adif 

/ 

calculate error 

4523 

/3772 


f mul 

i adif 

/ 

square it and 

4524 

/ 1356 


f add 

tot 

/ 

add to total 

4525 

/7356 


f put 

tot 



4526 

/00e0 


f ext 




4527 

/ 1361 


tad 

m4 

/ 

prepare to 

4530 

/3046 


dca 

ent 

/ 

type 

4531 

/ 1351 


tad 

am 7 

/ 

7 variables 

4532 

/3010 


dca 

1 0 



4533 

/ 1370 


tad 

m7 



4534 

/3047 


dca 

ch 

/ 

set counter 

4535 

/ 1410 


tad i 

. 10 

/ 

get next address 

4536 

/3034 


dca 

tm 

/ 

store it 

4537 

/4476 


f pnt 




4540 

/0434 


f get 

i tm 

/ 

get next variable 

4541 

/0000 


f ext 




4542 

/4746 


jms i 

. atypef / 

type it 

4543 

/2047 


isz 

ch 

/ 

check count 

4544 

/5335 


jmp 

.-7 



4545 

/5305 


jmp 

prt2+13 / 

return to get next dat 
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4546 

/7166 

atypef, 



4547 

/7300 

aw/xyz, 



4550 

/4600 

aply. 



4551 

/4570 

am 7^ 



4552 

/6200 

afout, 



4553 

/0077 

p77. 



4554 

/4753 

azro, 



4555 

/7550 

aone, 



4556 

/ 1736 

tot, 



4557 

/ 1403 




4560 

/4265 




4561 

/7774 

m4, 



4562 

/ 1755 

dif, 



4563 

/0276 




4564 

/2400 




4565 

/2056 

wc, 



4566 

/0120 




4567 

/5745 




4570 

/7771 

m7, 



4571 

/4565 

aw/c, 



4572 

/4562 

adif, 



4573 

/4467 

anum, 



4574 

/7564 

adw> 



4575 

/7567 

adx, 



4576 

/7572 

ady, 



4577 

/7575 

adz. 



4600 

/4512 

ply. 

Q 


4601 

/3044 


dca 

cy 

4602 

/ 1044 


tad 

cy 

4603 

/3010 


dca 

1 D 

4604 

/ 1344 


tad 

here 

4605 

/3ei 1 


dca 

11 

4606 

/ 1337 


tad 

m6 

4607 

/3043 


dca 

cx 

4610 

/ 1410 


tad i 10 

4611 

/341 1 


dca : 

L 11 

4612 

/2043 


isz 

cx 

4613 

/5210 


jmp 

• -3 

4614 

/7307 


7307 


4615 

/ 1044 


tad 

cy 

4616 

/ 1352 


tad 

n5 

4617 

/3341 


dca 

cn 

4620 

/ 1336 


tad 

ax 

4621 

/3340 


dca 

aout 

4622 

/7346 


7346 


4623 

/3351 


dca 

n4 

4624 

/ 1600 


tad 

i ply 

4625 

/3356 

bck, 

dca 

f X 

4626 

/4476 


f pnt 


4627 

/e756 


f get 

i f X 

4630 

/0000 


f ext 


4631 

/ 1351 


tad 

n4 

4632 

/ 1345 


tad 

nO 

4633 

/77 10 


spa 

cla 

4634 

/4743 


jms 

i alog 

4635 

/4476 


f pnt 


4636 

/7740 


f put 

i aout 

4637 

/0000 


f ext 



/ polynomial expansion routine 
/ to calculate a value for w 
/ that matches given values of 
/ x, y, &. z. Set index register 
/ get address where parameters 
/ are to be stored in this routine. 
/ set counter 

/ move parameters into this routine 
/ check counter 

/ set cn at address of the last 
/ regression coefficient. 

/ set aout 2 = to address where 
/ X is to be stored in this routine 
/ -3 

/ set counter 
/ get address of x from 
/ calling routine and store in fx. 

/ get next variable from calling 
/ program. 

/ check to see if we are 
/ to take the logarithm 
/ skip if no 
/ take the logarithm 

/ store the variable in this 
/ routine. 




2D 


4640 

/7325 


7325 


4641 

/ 1340 


tad 

aout 

4642 

/3340 


dca 

aout 

4643 

/7325 


7325 


4644 

/ 1356 


tad 

fx 

4645 

/2351 


isz 

n4 

4646 

/5225 


jmp 

bck 

4647 

/4476 


f pnt, 

, 

4650 

y"0353 


f get 

zro 

4651 

/7364 


f put 

f z 

4652 

70000 


f ext 


4653 

/ 1350 


tad 

n3 

4654 

73045 


dca 

cz 

4655 

7 1347 


tad 

n2 

4656 

73044 


dca 

cy 

4657 

74476 


f pnt 


4660 

70375 


f get 

'z 

4661 

73364 


fmul 

fz 

4662 

77364 


f put 

fz 

4663 

70353 


f get 

zro 

4664 

77361 


f put 

fy 

4665 

70000 


f ext 


4666 

7 1346 

iy. 

tad 

n1 

4667 

73043 


dca 

cx 

4670 

74476 


f pnt 


4671 

70372 


f get 

y 

4672 

73361 


fmul 

^y 

4673 

77361 


f put 

fy 

4674 

70353 


f get 

zro 

4675 

77356 


f put 

fx 

4676 

76600 


f skp 


4677 

74476 

lx, 

f pnt 


4700 

70367 


f get 

X 

4701 

73356 


fmul 

f X 

4702 

7 1741 


f add 

i cn 

4703 

77356 


f put 

fx 

47 04 

7 0000 


f ext 


4705 

77346 


7346 


4706 

7 1341 


tad 

cn 

4707 

73341 


dca 

cn 

47 10 

72043 


isz 

cx 

471 1 

75277 


jmp 

lx 

4712 

74476 


f pnt 


47 13 

70361 


fget 


4714 

7 1356 


f add 

f X 

4715 

77361 


f put 

^^y 

47 16 

70000 


f ext 


4717 

72044 


isz 

cy 

4720 

75266 


jmp 

ly 

4721 

74476 


f pnt 


4722 

70364 


fget 

fz 

4723 

7 1361 


f add 

fy 

4724 

77364 


f put 

fz 

4725 

70000 


f ext 


47 26 

72045 


isz 

cz 

4727 

75255 


jmp 

Iz 


/ increment aout by 3 

/ increment fx 
/ check count 

/ go back for next variable 


/ expand the polynomial 
/ equation. 



4730 

/7346 


7346 

4731 

/ 1345 


tad nO 

4732 

/7750 


spa sna cla 

4733 

/4742 


jms i aantl 

4734 

/2200 


isz ply 

4735 

/5600 


jmp i ply 

4736 

/4767 

ax, 


4737 

/7772 

m6, 


4740 

/5000 

aout, 


4741 

/0103 

cn, 


4742 

/5271 

aantl, 


4743 

/6661 

alog, 


4744 

/4744 

here, 


4745 

/0001 

nOi, 


4746 

/7775 

n1 , 


4747 

/7776 

n2, 


47 50 

/7777 

n3, 


4751 

/000e 

n4, 


4752 

/0022 

n5, 


4753 

70000 

zro, 


4754 

/0000 



4755 

/0000 



4756 

72024 

fx, 


4757 

77637 



4760 

70743 



4761 

72026 



4762 

72715 


4763 

73722 



4764 

72026 



4765 

72715 



4766 

73722 



4767 

72034 

X, 


4770 

73034 



4771 

76177 



4772 

72005 



4773 

74271 



4774 

7031 1 



4775 

72014 



4776 

70000 



4777 

70000 



0043/1744 

cx, 


0044 

70000 

cy, 


0045 

70000 

cz, 



/ check to see if we are to 
/ take the antilogarithm of 
/ the calculated value of w. 
/ yes 

/ increment return address 
/ return. 
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5000 

/0204 

mes, 

0 


5001 

/ 1217 


tad 

m1 

5002 

/3017 


dca 

1 7 

5003 

/ 1417 


tad 

i 17 

5004 

/3306 


dca 

chr 

5005 

/ 1306 


tad 

chr 

5006 

/2306 


isz 

chr 

5007 

/7410 


skp 


5010 

/5600 


jmp 

i mes 

5011 

/2306 


isz 

chr 

5012 

/7410 


skp 


5013 

/522G 


jmp 

. + 5 

5014 

/4616 


jms 

i apw 

5015 

/5203 


jmp 

mas+ 3 

5616 

/44@0 

apw, 



5017 

/7777 

ml , 



5020 

/4266 


jms 

cr 

5021 

/5203 


jfnp 

mes+3 

5043/6260 

afout, 



5044 

/4620 

s3, 



5045 

/4206 

s2, 



5046 

/0004 

s1 , 



5047 

/2427 

sO, 



5050 

/4453 

low. 



5051 

/4524 

high, 



5052 

/5405 

ret. 



5053 

/0215 

P215, 



5054 

/54e0 

afpnt, 



5055 

/63e5 

printc 

f 


5056 

/0212 

p212, 



5057 

/0240 

p240. 



5060 

/5073 

care, 

0 


5661 

/ 1253 


tad 

p21 5 

5062 

/4655 


jsm 

i prin 

5063 

/ 1256 


tad 

p212 

5064 

/4655 


jms 

i prin 

5065 

/5660 


jmp 

i care 

5066 

/4423 

cr, 

G 


5067 

/70G1 


iac 


5070 

/7650 


sna 

cla 

5071 

/5275 


jnip 


5072 

/4260 


jms 

care 

5073 

/2377 


isz 

Ic 

5074 

/ 5666 


jnip 

i cr 

5075 

/ 1377 


tad 

Ic 

5076 

/1375 


tad 

m6 

5077 

/3377 


dca 

Ic 

5106 

/4260 


jms 

care 

5101 

/2377 


isz 

Ic 

5102 

/5300 


jmp 

• -2 

5103 

/ 1376 


tad 

m74 

5104 

/3377 


dca 

Ic 

5105 

/5666 


jmp 

i cr 

5106 

/001 1 

chr, 



5167 

/7774 

m4, 




/ message subroutine 
/ subtract 1 from address 
/ set index register 
/ get next word 
/ save word 

/ is it eof ? 

/ yes. return, 

/ no. is it a carriage return ? 
/ yes. 

/ no. print word. 

/ go back for next word. 

/ address of print word. 

/ return carriage. 

/ go back for next word, 

/ address of output routine. 


/ low end of trace range. 

/ high end of trace range 
/ return address for trace. 

/ address of fpnt 
/ address of print character. 


/ basic carriage return 


/ carriage return 
/if -I » turn page. 


/ check line count. 

/ return if not empty 

/ count 6 lines at 
/ end of page 

/ count to start of next 
/ page. 

/ reset line count at 60. 
/ return 
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UTILITY PROGRAMS 


5110 

/5152 

poa, 

0 



/ 

print octal accumulator 

51 11 

77004 

ral 





5112 

73306 


dca 

chr 


/ 


5113 

71307' 


tad 

m4 


set count for four 

51 14 

73046 


dca 

cnt 


/ 

characters. 

51 15 

7 1306 


tad 

chr 




5116 

77004 


ral 


/ 

sh 

lift next character 

51 17 

77006 


rtl 



/ 

into position 

5120 

73306 


dica 

chr 


/ 


5121 

71306 


tad 

chr 


save character 

5122 

70332 


and 

P7 


/ 

mask digit 

5123 

71333 


tad 

p260 


/ 

make it printable 

5124 

74655 


jms 

i printc 


/ 

type this digit 

5125 

72046 


is 2' 

cnt 


/ 

check count 

5126 

75315 


jmp 

.-11 




5127 

71257 


tad 

p240 


/ 

type space 

5130 

74655 


jms 

i printc 


/ 

when done 

5131 

7 5710 


jmp 

i poa 


/ 

and return 

5132 

70007 

p7, 






5133 

70260 

p260, 




/ 


5134 

71654 

trace, 

tad 

i afpnt 


get address 

5135 

77041 


cia 



/ 

make it negative 

5136 

73247 


dca 

sO 


/ 

save 

5137 

7 1250 


tad 

low 




5140 

71247 


tad 

sO 




5141 

77740 


sma 

sza cia 




5142 

75373 


jmp 

dwn 


/ 

return if less than low 

5143 

71251 


tad 

high 




5144 

71247 


tad 

&Q 


/ 

return if beyond high 

5145 

77710 


spa 

cia 




5146 

75373 


jmp 

dwn 




5147 

74266 


jms 

cr 


/ 

carriage return 

5150 

71654 


tad 

i afpnt 

/ 

get address of instruction 

5151 

74310 


jms 

poa 


/ 

print octal address 

5152 

7 1053 


tad 

ac3 


/ 


5153 

73244 


dca 

s3 


sav fac 

5154 

7 1054 


tad 

ac2 




5155 

7 3245 


dca 

s2 




5156 

7 1055 


tad 

ad 




5157 

73246 


dca 

s1 




5160 

71061 


tad 

exp 




5161 

7 3247 


dca 

sO 




5162 

74643 


jms 

i fout 


/ 

type fac 

5163 

7 1247 


tad 

sO 


/ 

restore fac 

5164 

73061 


dca 

exp 




5165 

71246 


tad 

s1 




5166 

73055 


dca 

ad 




5167 

7 1245 


tad 

s2 




5170 

73054 


dca 

ac2 




5171 

71244 


tad 

s3 




5172 

73053 


dca 

ac3 




5173 

72654 

dwn, 

isz 

i afpnt 


/ 

increment addr. 

5174 

75652 


jmp 

i ret 


/ 

return 

5175 

77772 

m6, 






5176 

77704 

m74, 






5177 

77712 

Ic, 










POOR 


MAN‘5 


EDITOR 


This program will copy on the typewriter, any tape read by the 
keyboard. If you set a character code in the sw/itch register , 
the computer will halt before the selected character is typed , 


5200 

/6040 

pme, 

tfl 


5201 

/6032 


kcc 


5202 

/6031 


ksf 

> 

5203 

/5202 


jmp 

.-1 

5204 

/6034 


krs 


5205 

/3047 


dca 

ch 

5206 

/7604 


las 


5207 

/7710 


spa 

cla 

5210 

/5217 


jmp 

. + 7 

5211 

/7604 


las 


5212 

/0227 


and 

p777 

5213 

/7041 


cia 


5214 

/ 1047 


tad 

ch 

5215 

/7640 


sza 

cla 

5216 

/5222 


jmp 

.+4 

5217 

/7402 


hit 


5220 

/7410 


skp 


5221 

/5200 


jmp 

pme 

5222 

/ 1047 


tad 

ch 

5223 

/6041 


tsf 

.-1 

5224 

/5223 


jmp 

5225 

/6046 


tls 


5226 

/5201 


jmp 

pme+1 

5227 

/0777 

p777. 




If you wish to delete the selected character, depress EXAM , CONT. 

If you wish to select every character, set 5W0 = 1 

For example, if you want to copy a PAL source tape 
with comments to be added at end of each line: 

Put 5200 in address register; then set 5W 3-11 = 215. Insert 

tape in reader, and depress CONT , 

The computer will halt at end of each line. Turn reader off and 
depress EXAM , CONT. Type desired comment followed by a carriage 
return. Turn reader on. Repeat for each line. You must insert the 
carriage return, even when no comment is required; unless you depress 
CONT, but not EXAM. 
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5233 

/0203 p203, 





5234 

/5277 fix, 

0 

k 



5235 

/4664 

jms 

i aar1 



5236 

/4664 

jms 

i aar1 



5237 

/4664 

jms 

i aar1 



5240 

/ 1061 

tad 

exp 



5241 

/7041 

cia 




5242 

/ 1355 

tad 

p200 



5243 

/7500 

sma 




5244 

/5251 

jmp 

• + 5 



5245 

/3046 

dca 

cnt 

/ 

no of left shifts required 

5246 

/4665 

jms 

i a all 

/ 

to justify decimal point. 

5247 

/2046 

isz 

cnt 



5250 

/5246 

jmp 

.-2 



5251 

/7600 

cla 




5252 

/3054 

dca 

ac2 

/ 

zero low order 

5253 

/3053 

dca 

ac3 

/ 

accumulator 

5254 

/ 1233 

taci 

p203 



5255 

/3061 

dca 

exp 

/ 

exponent for fixed part. 

52 56 

/1051 

tad 

sgn 



5257 

/7104 

cll 

ral 

/ 

put sign in link 

5260 

/ 1055 

tad 

ac1 



5261 

/7430 

szl 




5262 

/7041 

cia 


/ 

accumulator now holds fixed 

5263 

/5634 

jmp 

i fix 

/ 

part,, return. 

5264 

/ 5600 aar1 , 





5265 

/6050 aal1, 





5266 

/0201 p20l. 





5267 

/7777m1 , 





5270 

/5400 afpnt, 
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(X ) 

anti = ANTI LOG„ = e-.^ ' 

6 


5271 

/ 1664 

anti, 

0 


5272 

/4670 


jms i 

afpnt 

5273 

/3375 


f mul 

log2e 

5274 

/7040 


f put 

f 

5275 

/0000 


fext 


5276 

/4234 


jms 

fix 

5277 

/7510 


spa: 


5300 

/ 1267 


tad 

m1 

5301 

/3047 


dca 

ch 

5302 

/4670 


jms i 

afpnt 

5303 

/2040 


f sub 

f 

5304 

/7040 


f put 

f 

5305 

/6740 


f sgt 


5306 

/5312 


f jmp 

,+4 

5307 

/0372 


fget 

one 

5310 

/2040 


f sub 

f 

531 1 

/7040 


f put 

f 

5312 

/0000 


f ext 


5313 

/3051 


dca 

sgn 

5314 

/ 1047 


tad 

ch 

5315 

/ 1266 


tad 

p20l 

5316 

/3061 


dca 

exp 

5317 

/ 1354 


tad 

p4 

5320 

/3055 


dca 

ac1 

5321 

/3054 


dca 

ac2 

5322 

/3053 


dca 

ac3 

5323 

/4670 


jms i 

. afpnt 

5324 

/7043 


f put 

int 

5325 

/0040 


f get 

f 

5326 

/0000 


fext 


5327 

/3051 


dca 

sgn 

5330 

/4670 


jms i 

. afpnt 

5331 

/3356 


f mul 

ln2D2 

5332 

/7035 


f put 

y 

5333 

/3035 


f mul 

y 

5334 

/ 1361 


tad 

b1 

5335 

77040 


f put 

f 

5336 

70364 


fget 

a 

5337 

74040 


f div 

f 

5340 

72035 


fsub 

y 

5341 

7 1367 


f add 

aO 

5342 

77040 


f put 

f 

5343 

70035 


fget 

y 

5344 

7 1035 


f add 

y 

5345 

74040 


f div 

f 

5346 

7 1372 


fadd 

one 

5347 

77040 


f put 

f 

5350 

73040 


f mul 

f 

5351 

73043 


f mul 

int 

5352 

70000 


fext 


5353 

75671 


jmp i 

. anti 

5354 

70004 

p4, 



5355 

70200 

p2Q0, 




/ store characteristic 
/ fix - f s= 

/ negative mantissa 

/ compliment if diff. is 
/ positive, (mantissa is neg) 


/ make mantissa positive. 


/ fac now. contains 
/ 2 to power 


/ using formulas on 
/ page 8-26 .of 

/ ’^introduction to prog.”* 
/ we now calculate the 
/ r-th power of 2, 



27 


5356 /1775 

5357 /4271 

5360 /0273 

5361 /2067 

5362 /4056 

5363 /1323 

5364 /6124 

5365 /5471 

5366 /5707 

5367 /2046 

5370 /0036 

5371 /6017 

5372 /2014 

5373 /0000 

5374 /0000 

5375 /2015 

5376 /6125 

5377 /0722 


ln2o2, 

/ .34657359 

b1 , 

/ 61 O.O 9 OI 9 I 

a1 , 

/-601.80427 

aD 

/ 12.015017 

one. 

/ 1 . 

log2e 

/ 1 .442695 
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6502 

/0262 

rint f 

0 


6503 

/3061 


dca 

exp 

6504 

/3062 


dca 

dot 

6505 

/3051 


dca 

sgn 

6506 

/3053 


dca 

ac3 

6507 

/3054 


dca 

ac2 

6510 

/3055 


dca 

a cl 

6511 

/7240 


cla 

cma 

6512 

/3046 


dca 

cnt 

6513 

/6031 

loop, 

ksf 


6514 

/5313 


jrnp 

.-1 

6515 

/6036 


krb 


6516 

/3047 


dca 

ch 

6517 

/1047 


tad 

ch 

6520 

/1357 


tad 

m2l 2 

6521 

/5610 


jmp 

i a74l0 

6522 

/1047 


tad 

ch 

6523 

/1363 


tad 

m256 

6524 

/7440 


sza 


6525 

/5331 


jmp 

.+4 

6526 

/7240 


cla 

cma 

6527 

/3062 


dca 

dot 

6530 

/5313 


jmp 

loop 

6531 

/7001 


iac 


6532 

/7440 


sza 


6533 

/5337 


jmp 

.+4 

6534 

/7330 


7330 

6535 

/3051 


dca 

sgn 

6536 

/5313 


jmp 

loop 

6537 

/1365 


tad 

m3 

6540 

/7510 


spa 


6541 

/5313 


jmp 

loop 

6542 

/1364 


tad 

ml 2 

6543 

/7700 


sma 

cla 

6544 

/5702 


jmp 

i rint 

6545 

/4760 


jms 

i ami 0 

6546 

/1047 


tad 

ch 

6547 

/ 1362 


tad 

m260 

6550 

/3056 


dca 

op3 

6551 

/3057 


dca 

op2 

6552 

/3060 


dca 

opi 

6553 

/4761 


jms 

i aoadd 

6554 

/1062 


tad 

dot 

6555 

/1046 


tad 

cnt 

6556 

/5312 


jmp 

loop-1 

6557 

/7566 




6560 

/6153 

am1 0, 



6561 

/5633 

aoadd, 



6562 

/7520 

m26D, 



6563 

/7522 

m256, 



6564 

/7766 

m1 2, 



6565 

/7775 

m3, 




/ read integer 
/ initialize variables 


/ is it a linefeed ? 

/ patch to count line feeds 

/ decimal point ? 

/no, 

/ yes. 

/ dot = -1 


/ minus sign ? 

/ no. 

/ 4000 

/ set sign = 4000 

/ greater than 260 ? 

/ no, get next digit. 

/ yes 

/ less than ?72 ? 

/ return if no. 

/ yes. mult, fac by 10, 
/ recover character. 

/ make it a digit 
/ put in operand. 


/ fac = fac + operand 
/ decriment count 

/ if a point (.) was encountered 
/ get next digit. 

/ minus 21 2 



f 


« * 
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6600 

/0033 

fin, 

□ 


/ 

6601 

/4657 


jms 

i arint 

/ 

6602 

/1255 


tad 

p233 

/ 

6603 

/3061 


dca 

exp 

/ 

6604 

/1055 


tad 

ac1 

6605 

/1254 


tad 

m4 

/ 

6606 

/7700 


sma 

cla 

/ 

6607 

/5214 


jmp 

dwn 

6610 

/4650 


jms 

i aall 


6611 

/7240 


cla 

cm a 


6612 

/1061 


tad 

exp 

/ 

6613 

/5203 


jmp 

.-10 

6614 

/4342 

dwn. 

jms 

mar 

/ 

6615 

/2046 


isz 

cnt 

/ 

6616 

/7410 


skp 


/ 

6617 

/5600 


3^9 

i fin 

/ 

6620 

/1253 


tad 

ml 1 


6621 

/3261 


dca 

cl 1 

/ 

6622 

/1053 


tad 

ac3 

/ 

6623 

/3056 


dca 

□ Pp3 


6624 

/3053 


dca 

ac3 


6625 

/1054 


tad 

ac2 


6626 

/3057 


dca 

op?2 


6627 

/3054 


dca 

ac2 


6630 

/1055 


tad 

ac1 


6631 

/3060 


dca 

opi 


6632 

/3055 


dca 

ac1 


6633 

/4651 

bkk, 

jms 

i aarl 

/ 

6634 

/4651 


jms 

i aarl 


6635 

/4651 


jms 

i aarl 


6636 

/4652 


jms 

i aoadd 


6637 

/4651 


jms 

i aarl 


6640 

/4652 


jms 

i aoadd 


6641 

/2261 


isz 

cl 1 

/ 

6642 

/5233 


jmp 

bck 

6643 

/1254 


tad 

m4 

/ 

6644 

/1061 


tad 

exp 

/ 

6645 

/3061 


dca 

exp 


6646 

/5214 


jmp 

dwn 


6647 

/0000 





6650 

/6050 

aal1 , 




6651 

/5600 

aar1 , 




6652 

/5633 

aoaddl, 




6653 

/7767 

mil 1 , 




6654 

/7774 

m4, 




6655 

/0233 

p233, 




6656 

/0000 





6657 

/6502 

arint 




6660 

/7770 

Tn1 Q 





floating input routine 
read integer 
initialize exponent 

shift fac left 1 to fill 
low order dig. in ac1 

when finished 


loob back 

mo^e fac right to put only 

one digit in ad . 

count digits right of dec, point 

return when done. 

c11 =11 

move fac to operand and clear 


miultiply fac by 16/10 


loop bck 9 times, 
subtract 4 from exponent 
each time 





6661 

/0000 log, 

0 


6662 

/4752 

jms i 

afpnt 

6663 

/7040 

f put 

f 

6664 

/0000 

f ext 


6665 

/ 1061 

tad 

exp 

6666 

/7550 

spa sna 

6667 

/7402 

hit 


6670 

/1353 

tad 

mZDO 

6671 

/3200 

dca 

tm 

6672 

/1354 

tad 

p2D3 

6673 

/3061 

dca 

exp 

6674 

/7330 

g4000 

6675 

/0200 

and 

tm 

6676 

/3051 

dca 

sgn 

6677 

/3053 

dca 

3c3 

6700 

/3054 

dca 3 c2 

6701 

/1200 

tad 

tm 

6702 

/7510 

spa 


6703 

/7041 

cla 


6704 

/3055 

dca 

ac1 

6705 

/4342 

jms 

mar 

6706 

/4752 

jms i 

afpnt 

6707 

/7043 

f put 

int 

6710 

/0040 

fget 

f 

6711 

/0000 

f ext 


6712 

/1355 

tad 

p200 

6713 

/3061 

dca 

exp 

6714 

/3051 

dca 

sgn 

6715 

/4752 

jms i 

afpnt 

6716 

/7040 

f put 

f 

6717 

/ 1364 

f add 

sqh 

6720 

/7035 

f put 

z 

6721 

/0040 

fget 

f 

6722 

/2364 

f sub 

sqh 

6723 

/4035 

f div 

z 

6724 

/7035 

f put 

z 

6725 

/3035 

f mul 

z 

6726 

/7040 

f put 

f 

6727 

/0375 

f get 

c5 

6730 

/3040 

f mul 

f 

6731 

/ 1372 

fadd 

c3 

6732 

/3040 

f mul 

f 

6733 

/ 1367 

fadd 

cl 

6734 

/3035 

f mul 

z 

6735 

/2356 

f sub 

phlf 

6736 

/ 1043 

fadd 

int 

6737 

/3361 

f mul 

ln2 

6740 

/0000 

f ext 


6741 

/5661 

jmp i 

. log 


/ calculate naperian logarithm 

/ save argument 

/ get exponent of arg. 

/ halt if zero or neg. 

/ change exp to + or - . 

/ initialize exponent 
/ of integer part 

/ insert proper sign 

/ put original exponent 
/ in fac 

/ justify fac right 
/ store integer part 

/ mantissa exponent is zero. 

/ sign is also zero 

/ store mantissa 

/ z = f + sqrt(.5) 

/ see "'Introduction to Prog," 

/ page 8-27 . 

/ z = (f - sqrt( .5) )/ z 

/ f = z**2 

/ fac = ((c5*f + c3)*f + 

/ c1)*z - .5 + int 

/ convert to naperian base 

/ return 



6742 

/6615 

mar, 0 


6743 

/1055 

tad 

ac1 

6744 

/ 1260 

tad 

ml 0 

6745 

/7710 

spa 

cla 

6746 

/5742 

jn»P 

i mar 

6747 

74652 

jms 

i aarl 

6750 

72061 

isz 

exp 

6751 

75343 

jmp 

.-6 

6752 

75400 

afpnt, 


6753 

77600 

m200. 


6754 

70203 

p203, 


6755 

70200 

p200, 


6756 

7 1777 

phlf, 


6757 

77777 



6760 

77774 

m4, 


6761 

72005 

ln2, 


6762 

74271 



6763 

70276 



6764 

72005 

s>qh, 


6765 

75202 



6766 

73627 



6767 

72025 

cl , 


6770 

76125 



6771 

70776 



6772 

72007 

c3, 


6773 

7 5421 



6774 

7 3571 



6775 

72004 

c5, 


6776 

76253 



6777 

72516 




/ / move fac right to empty 
/ all bits left of bit 9 in ac1 

/ are we through ? 

/ yes, return 
/ no. shift fac right. 

/ correct ecponent 

/ address of f.p,inter. 

/ .5 

/ .6931472 
/ .7071 068 
/ 2.885391 3 
/ .9614706 
/ .59897865 
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SOLVE SIMULTANEOUS EQUATIONS 


■fortran program 


0073 

/5110 

po9» 


do 10 k=1,n 

0074 

/5067 

cr. 


f = X (k, k ) 

0075 

/6305 

type. 


do 10 j = k,n 

0076 

/5400 

f ppt. 


x{k,j+1) + x(k,j+t)/ f 

0077 

/5134 

link, 


do 10 i = 1 , n 

0100 

0101 

^7775 

/0000 

n. 


if( i ,eq. k) go to 10 
x{i,j+1) = x(i,j+1) - ; 

0102 

0103 

0104 

0105 

/0000 
/ 0000 
/0000 
/001 1 

m, 

1 0 

cointinue 


7000 

/ 0000 

ic. 




7001 

/0000 

jc. 




7002 

/0302 

kj » 




7003 

/0274 

ij. 




7004 

/0000 

kc. 




7005 

/0252 

ok, 




7006 

/0260 

kk , 




7007 

/0263 

ik. 




7010 

/0000 

/ 




‘7011 

/0236 

at, 




7012 

/5023 

sse, 

□ 

/ 

solve simultaneous 

7013 

/1211 


tad 

at 


7014 

/3206 


dca 

kk / 

initialize addresses 

7015 

/1100 


tadi 

n 


7016 

/3204 


dca 

kc 


7017 

/1211 


tad 

at 


7020 

/3205 

k, 

dca 

ok 


7021 

/1205 


tad 

ok 


7022 

/ 1105 


tad 

m 


7023 

/3203 


dca 

ij 


7024 

/4476 


jms 

i ifpnt 


7025 

/0606 


get 

i kk / 

get term on diagonal 

7026 

/7040 


f put 

f 


7027 

/0000 


f BXt 



7030 

/1206 


tad 

kk 


7031 

/1105 


tad 

m 


7032 

/3202 


dca 

kj 


7033 

7 1204 


tad 

kc 


7034 

/3201 


dca 

jc 


7035 

/4476 

j » 

jms 

i ifpnt 


7036 

70602 


f get 

i kj 


7037 

74040 


fdiv 

f / divide key row by term on 

7040 

77602 


f put 

i kj 


7041 

7 0000 


T 



7042 

71205 


tad 

ok 


7043 

73207 


dca 

ik 


7044 

71100 


tad 

n 


7045 

7 3200 


dca 

ic 



(k,j+1)* x{i,k) 


diagonal 



7046 

/ 1200 

« • 

tad 

ic 

7047 

/7041 


cia 


7050 

/ 1204 


tad 

kc 

7051 

/7650 


jmp 

cla 

7052 

/5263 


jmp 

sk 

7053 

/4476 


f pnt 


7054 

/0602 


fget 

i kj 

7055 

/3607 


f mul 

i ik 

7056 

/7043 


f put 

int 

7057 

/0603 


fget 

i ij 

7060 

72043 


f sub 

int 

7061 

77603 


f put 

i ij 

7062 

70000 


f ext 


7063 

77325 

sk, 

7325 


7064 

7 1203 


tad 

ij 

7065 

73203 


dca 

ij 

7066 

77325 


7325 


7067 

7 1207 


tad 

ik 

7070 

73207 


dca 

ik 

7071 

72200 


isz 

ic 

7072 

75246 


jniP 

ii 

7073 

7 1105 


tad 

m 

7074 

7 1202 


tad 

kj 

7075 

73202 


dca 

kj 

7076 

72201 


isz 

jc 

7077 

75235 


jmp 

j 

7100 

72204 


isz 

kc 

7101 

77410 


skp 


7102 

75315 


jmp 

ma 

7103 

77325 


7325 


7104 

7 1105 


tad 

m 

7105 

7 1206 


tad 

kk 

7106 

7 3206 


dca 

kk 

7107 

7 1205 


tad 

ok 

7110 

7 1105 


tad 

m 

71 1 1 

75220 


jmp 

k 

71 12 

77772 

m6, 



7113 

7 0000 




71 14 

7 0000 





/ skip operation if i = k 

/ generate zeros on key column 


/ 3 

/ increment ij 
/ 3 

/ increment ik 
/ iterate down column 

/ increment kj 
/ iterate across row 

/ 3 

/ increment kk 

/ go back to sweep matrix again 
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7U5 

/1105 

ma, 

tad 

m / move 

7116 

/7041 


cla 


7117 

/3204 


dca 

kc 

7120 

/7305 


7305 

1 1 -' / 2 

7121 

/1206 


tad 

kk 

7122 

/3010 


dca 

10 

7123 

/1364 


tad 

pi 05 

7124 

/3011 


dca 

11 

7125 

/1410 


tad 

i 10 

7126 

/3411 


dca 

i 11 

7127 

/2204 


isz 

kc 

7130 

/5325 



.-3 

7131 

/7600 


cla 


7132 

/ 1312. 


tad 

mi6 

7133 

>'3204 


dca 

kc 

7134 

>'4474 


jms 

i cr 

7135 

/1363 


tad 

p77 

7136 

/3010 


dca 

10 

7137 

/1410 


tad 

i 10 

7140 

/4473 


jmji 

poa / print 

7141 

/2204 


is 2 

kc 

7142 

75337 


jnip 

.-3 

7143 

77240 


cla 

cma 

7144 

73046 


dca 

cnt 

7145 

71100 


tad 

n 

7146 

73204 


dca 

kc 

7147 

72010 


isz 

1 0 

7150 

74476 


jms i fpnt 

7151 

70410 


fget 

i 10/ print 

7152 

70000 


f ext 


7153 

72010 


isz 

1 0 

7154 

72010 


isz 

1 0 

7155 

74366 


jms 

i typef 

7156 

72204 


isz 

kc 

7157 

75347 


jmp 

.-10 

7160 

75612 


jmp 

i sse 

7161 

76200 

af out 

f . 


7162 

77774 

m4, 



7163 

70077 

P77, 



7164 

70105 

pl05. 



7165 

70240 

p240, 



7166 

77156 

typef 

,0 

/ type 

7167 

72046 


isz 

cnt 

7170 

75374 


jmp 

.+4 

7171 

74474 


jms 

i cr 

7172 

71362 


tad 

m4 

7173 

73046 


dca 

cnt 

7174 

74761 


jms 

af out 

7175 

71365 


tad 

p240 

7176 

74475 


jms 

i type 

7177 

75766 


jmp 

i typef 


1 06 


/ print six parameters in octal 
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TRACE OF 55E Solve simultaneous equations 

START 5020G 


7061 

1 .000000£-t-00 

7025 

2.000000E+00 

7025 

1.500000E+00 

7026 

2. 000000E+00 

7026 

1.500000E-*-00 

7036 

3.000000E+00 

7036 

1.000000E+00 

7037 

1.500000E+00 

7037 

6.666667E-01 

7040 

1.500000E+00 

7040 

6.666667E-01 

7054 

1.500000E+00 

7054 

6.666667E-01 

7055 

4.500000E-f00 

7055 

1.000000E+00 

7056 

4.500000E+00 

7056 

1.000000E+00 

7057 

6 . 000000E+00 

7057 

2.000000E+00 

7060 

1.500000E+00 

7060 

1.000000E-t-00 

7061 

1.500000E+00 

7061 

1.000000E-I-00 

7054 

1 . 500000E+00 

7054 

6.666667E-01 

7055 

6 * 000000E't'00 

7055 

6.666667E-01 

7056 

6 .000000E+00 

7056 

6.666667E-01 

7057 

7.000000E+00 

7057 

1.000000E+00 

7060 

1.000000E4-00 

7060 

3.333333E-01 

7061 

1.000000E+00 

7061 

3.333333E-01 

7036 

4. 000000E+00 

7036 

5. 000000E-01 

7037 

2, 000000E-f'00 

7037 

3. 333333E-01 

7040 

2 . 000000E-I-00 

7040 

3.333333E-01 

7054 

2. 000000E+00 

7054 

3. 333333E-01 

7055 

6 . 000000E-I-00 

7055 

5.000000E-01 

7056 

6 * 000000E‘<‘00 

7056 

5. 000000E-01 

7057 

7. 000000E+00 

7057 

5.000000E-01 

7060 

1 .000000E-t-00 

7060 

0 

7061 

1.000000E+00 

7061 

0 

7054 

2.000000E+00 

7054 

3.333333E-01 

7055 

8 .000000E+00 

7055 

3.333333E-01 

7056 

6 .000000E+00 

7056 

3. 333333E-01 

7057 

9* 000000E+00 

7057 

1.000000E+00 

7060 

1.000000E+00 

7060 

6*666667E-01 

7061 

1 .000000E-f00 

7061 

6.666667E-01 

7036 

1 .000000E+00 

7025 

3. 333333E-01 

7037 

5.000000E~ 

01 

7026 

3.333333E-01 

7040 

5.000000E- 

01 

7036 

6.666667E-01 

7054 

5.000000E- 

01 

7037 

2. 000000E4-00 

7055 

!•500000E+00 

7040 

2.000000E+00 

7056 

1•500000E+ 00 

7054 

2. 000000E-t-00 

7057 

2# 000000E’^00 

7055 

2* 000000E+00 

7060 

5.000000E- 

01 

7056 

2. 000000E+00 

7061 

5.000000E- 

01 

7057 

0 

7054 

5.000000E- 

01 

7060 

>2. 000000E+00 

7055 

2. 000000E-I-00 

7061 

-2.000000E+00 

7056 

2.000000E+00 

7054 

2.000000E+00 

7057 

3.000000E+00 

7055 

1.333333E+00 

7060 

1.000000E+00 

7056 

1.333333E+00 




7057 

3.333333E-01 




7060 

> 1.000000E+00 




7061 

- 1.000000E-*-00 



-2 

7775 

0000 0000 0000 

0000 0011 



7151 

-2.000000E+00- 

2.000000E+00 

2 

7151 

- 1.000000E+00- 

1.000000E+00 

3 

7151 

2.000000E+00 

2.000000E-I-00 



TEST MATRIX 


-1 
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MULTIPLE REGRESSION PROGRAM 


7200 

/6032 ra1 p 

kcc 


/ 

7201 

/6040 

tfl 



7202 

/4666 

jms 

i arint 

/ 

7203 

/1047 

tad 

ch 

/ 

7204 

/1373 

tad 

m305 

/ 

7205 

/7450 

sna 



7206 

/5222 

jmp 

j-1 A 

. » -• —r 

/ 

/ 

7207 

/7500 

sma 



7210 

/5202 

jmp 

ra1 +2 

/ 

7211 

/ 1372 

tad 

p5 

/ 

7212 

/7550 

spa 

sna 


7213 

75202 

jmp 

ra1 +2 

/ 

7214 

71371 

tad 

p31 00 

/ 

7215 

73220 

dca 

. + 3 


7216 

7 1053 

tad 

ac3 

/ 

7217 

77041 

cia 


/ 

7220 

73102 

3101 


/ 

7221 

7 5202 

jmp 

ra1 +2 

/ 

7222 

7 1 101 

tad 

n1 


7223 

77500 

sma 


/ 

7224 

77402 

hit 


/ 

7225 

77041 

cia 


/ 

7226 

73105 

dca 

nS 

/ 

7227 

71102 

tad 

n2 

/ 

7230 

74267 

jms 

mli 

/ 

7231 

7 1103 

tad 

n3 

/ 

7232 

74267 

jms 

mli 

/ 

7233 

7 1105 

tad 

n5 


7234 

77041 

cia 



7235 

73100 

dca 

nO 

/ 

7236 

77346 

7346 


/ 

7237 

74267 

jms 

mli 

/ 

7240 

7 1100 

tad 

nO 

/ 

7241 

7 1370 

tad 

p34 


7242 

77710 

spa 

cia 


7243 

77402 

hit 


/ 

7244 

71100 

tad 

nO 


7245 

7 3046 

dca 

cnt 

/ 

7246 

7 1105 

tad 

n5 

/ 

7247 

72046 

isz 

cnt 


7250 

7 5246 

jmp 

.-2 


7251 

71105 

tad 

n5 


7252 

77041 

cia 



7253 

73046 

dca 

cnt 

/ 

7254 

7 1367 

tad 

p235 

/ 

7255 

73010 

dca 

10 


7256 

73410 

dca 

i 10 

/ 

7257 

72046 

isz 

cnt 


7260 

75256 

jmp 

.-2 


7261 

74766 

jms 

i areg 

/ 

7262 

77300 

cll 

cia 


7263 

74765 

jms 

i asse 

/ 

7264 

77402 

hit 


/ 

7265 

75764 

jmpi 

. aprt2 

/ 


clear &. set flags 

get number 
test termination 
is it e? 

yes. end of parameters 

go back if greater than e. 
is it less than a. 

yes. go back for next char, 
no. fix up a store instruction 

get the integer 
make it negative 
store it. 

get next parameter. 

get first parameter, 
halt if positive 
make it positive 
load in multiplicand reg. 
get second parameter, 
multiply n1*n2 
get third parameter 
multiply it also 

-(n1»n2*n3) 

-3 

multiply prod by 3 
no of degrees of freedom 

halt if prod. .gt. 20 
set count = nO 

multiply by words per column 


count = no. of words in 
matrix, set index register 

clear entire matrix 


form regression matrix 

solve simultaneous equations 
halt when done 
jump to part 2 
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7266/6502 

arint, 



/ 

7267 

/7240 

mli, 

0 


/ 

7270 

/7500 


sma 



7271 

/7402 


hit 


/ 

7272 

/3046 


dca 

ent 


7273 

/1105 


tad 

n5 


7274 

/2046 


isz 

ent 


7275 

/5273 


jmp 

.-2 


7276 

/3105 


dca 

n5 

/ 

7277 

/5667 


jmp i 

mli 

/ 

7300 

/4506 

wxyz. 

0 


/ 

7301 

/4761 


jms i 

af in 

/ 

7302 

/1047 


tad 

ch 

/ 

7303 

/1363 


tad 

m333 


7304 

/7450 


sna 



7305 

/5326 


jmp i 

wxyz 

/ 

7306 

/7500 


sma 



7307 

/5700 


jmp i 

wxyz' 

7‘ 

7310 

/1372 


tad 

p5 

/ 

7311 

77550 


spa 

sna 

/ 

7312 

/5301 


jmp 

wxyz+l 

/ 

7313 

/1362 


tad 

p7773 

/ir 

7314 

73323 


dca 

ins 

/ 

7315 

77307 


7307 


/ 

7316 

71104 


tad 

n4 

/ 

7317 

71323 


tad 

ins 

/ 

7320 

77510 


spa 



7321 

74760 


jms i 

log 


7322 

74476 


f pnt 



7323 

77775 

ins, 

7775 


/ 

7324 

70000 


f ext 



7325 

75301 


jmp 

wxyz+1 

/ 

7326 

77000 


nop 


/ 

7327 

77410 


skp 


/ 

7330 

75301 


jmp 

wxy2+1 

/ 

7331 

72300 


isz 

wxyz 

/ 

7332 

75700 


jmp i 

. wxyz 

/ 

7333 

70000 





7334 

70000 





7335 

70000 





7336 

70000 





7337 

70000 





7340 

75177 

ale , 



/ 

7341 

77550 

aop e , 



/ 

7342 

74476 

ra. 

■fpnt 


/ 

7343 

70741 


f get 

i aone / 

7344 

77775 


f put 

i ax 

/ 

7345 

77776 


•fput 

i ay 

/ 

7346 

77777 


f put 

i az 


7347 

70000 


fext 



7350 

7 7240 


cl a cm a 


7351 

73740 


dca i ale 

/ 

7352 

77240 


cl a 

cma 

/ 

7353 

73102 


dca 

n2 


7354 

77240 


cl a 

cma 


7355 

73103 


dca 

n3 


7356 

73104 


dca 

n4 


7357 

75200 


jmp 

ra1 



address of read integer routine 
multiply integer routine 

halt if count pos. 


put product in parameter 5 
return 

read and store w, x, y, or z 
get floating point numb, from 
keyboard, test terminating charac 

end line if = 333 

eThci' of file i'f'"gf^eaYeY" tharr'TJsr' 

if = or less than 326, go back 

for next character, 

fix up a floating point point 

st. to store numb, in address 

specified by term, character 

4 

check n4 to determine if we 
are to take log. 


store number 

go back for next number, 

change tp a hit if you want to ed 

the input, EXAM, CONT will 

overwrite the current line, 

not end of file; so skip, 

and return. 


address of line count, 

address of one 

regress analysis program, 

set parameters at default 
values. 


set linecount at -1 to break 
page at start of problem. 
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7360 

/6661 

alog. 


/ 

address of log routine 

7361 

/6600 

af in, 


/ 

aHHrass of floating input rout. 

7362 

/7773 

p7773, 


/ 

plus 7773 

7363 

/7445 

m333, 


/ 

minus 333 

7364 

/4472 

aprt2. 


/ 

address of part two 

7365 

/7012 

a s s e, 


/ 

address of solve simul. equations 

7366 

/7413 

areg, 


/ 

address to build regression matri 

7367 

/0235 

p235, 


/ 

plus 235 

7370 

/0034 

p34, 


/ 

plus 34 

7371 

/3100 

p3l00, 



7372 

/0005 

p5, 




7373 

/7473 

m305, 




7374 

/7564 

aw, 


/ 

address of w 

7375 

/7567 

ax, 


/ 

address of x. 

7376 

/7572 

ay. 


/ 

address of y 

7377 

/7575 

az, 


/ 

address of z 

7400 

/0000 





7401 

/0000 





7402 

/6522 

a6522. 




7403 

/4474 


cr 



7404 

/5602 


jmp i a6522 


7405 

/1047 


tad ch 



7406 

/4475 


type 

/ 

this patch was 

7407 

/5602 


jmp i a6522 

7410 

/7640 


sza cia 

/ 

called from 6521 to 

7411 

/5205 


jmp .-4 

/ 

count line feeds entered 

7412 

/5203 


jmp .-6 

/ 

with the input. 

7413 

/7262 

reg, 

0 

/ 

routine to build regression matri: 

7414 

/4476 


f pnt 

/ 

enter floating point interpreter 

7415 

/0350 


fget one 


7416 

/7361 


fput VC 



7417 

/0000 


f ext 

/ 

exit floating point interpreter. 

7420 

/ 1341 


tad p106 


7421 

/3005 


dca uk 

/ 

address of subscripted var. u(k). 

7422 

/4743 


jms 1 awx\/z 

/ get w, X, y, z, from keyboard. 

7423 

/5613 


jmp 1 reg 


/ return if end of file 

7424 

/1103 


tad n3 


/ set up first row of matrix in i 

7425 

/3347 


dca cc 


/ set c count. 



7426 

/ 1102 

Ic, 

tad 

n2 

/ 

loop c 

7427 

/3346 


uta 

bu 

/ 

/ 

sat b count 39 

7430 

/4476 


f pnt 




7431 

/0350 


f get 

one 



7432 

/7356 


f put 

vb 



7433 

/6600 


f skp 




7434 

/4476 

lb. 

f pnt 


/ 

loop b 

7435 

/0350 


f get 

one 



7436 

/'7353 


f put 

va 



7437 

/0000 


f ext 




7440 

/1101 


tad 

n1 



7441 

/3345 


dca 

ac 

/ 

set a count 

7442 

/4476 


f pnt 




7443 

/0353 


f get 

va 



7444 

/3356 


fmul 

vb 



7445 

/3361 


fmul 

VC 


place value in u(k) 

7446 

/7405 


f put 

i uk 

/ 

7447 

/0353 


f get 

V3 



7450 

/3367 


fmul 

X 



7451 

/7353 


f put 

va 



7452 

/0000 


f B'Xt 




7453 

/7325 


7325 


/ 

3 

7454 

/1005 


tad 

uk 

/ 

increment to next tripple precision 

7455 

/3005 


dca 

uk 

/ 

word. 

7456 

/2345 


isz 

ac 

/ 

check a count 

7457 

/5242 


jn>P 

la 

/ 

repeat loop a 

7460 

/4476 


f pnt 




7461 

/0356 


fget 

vb 



7462 

/3372 


fmul 

y 



7463 

/7356 


f put 

vb 



7464 

/0000 


f ext 




7465 

/2346 


isz 

be 

/ 

check b count 

7466 

/5234 


jmp 

lb 

/ 

repeat loop b 

7467 

/4476 


f pnt 




7470 

/0361 


fget 

VC 



7471 

/3375 


fmul 

z 



7472 

/7361 


f put 

VC 



7473 

/0000 


f ext 




7474 

/2347 


isz 

cc 

/ 

check c count 

7475 

/5226 


jmp 

Ic 

/ 

repeat loop c 

7476 

/4476 


f pnt 




7477 

/0364 


fget 

u 



7500 

/7405 


f put 

i uk 

/ 

put w in last term of u(k) 

7501 

/0000 


f ext 

p236 



7502 

/1342 


tad 

/ 

set first address of matrix 

7503 

/3006 


dca 

xl 

/ 

in xl 

7504 

/7040 


cm a 


/ 

subtract 1 from negative count 

7505 

/I 100 


tad 

nO 

/ 

in npQ to process (n+1) terms. 

7506 

/3346 


dca 

me 



7507 

/ 1341 


tad 

p1 06 

/ 

set up first address of variable, u. 
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7510 

/3003 

ol, 

dca um 

/ store next address of u(m) 

7511 

/ 1100 


tad nO 

/ in um. to begin outer loop. 

7512 

/3345 


dca Ic 

/ set 1 count 

7513 

/ 1341 


tad p106 

/ initialize address 

75U 

/3007 


dca ul 

/ ul 

7515 

/4476 

il» 

f pnt 

/ inner loop 

7516 

/0407 


fget i ul 

/ get u(l) 

7517 

/3403 


fmul i urn 

/ multiply by u(m) 

7520 

/1406 


fadd i xl 

/ add to x(l) 

7521 

/7406 


fput i xl 

/ and update x(l) 

7522 

/0000 


f ext 

/ 3 

7523 

/7325 


7325 

7524 

/1007 


tad ul 

/ increment addiress ul 

7525 

/3007 


dca ul 

/ 3 

7526 

/7325 


7325 

7527 

/1006 


tad xl 

/ increment address xl 

7530 

/3006 


dca xl 


7531 

/2345 


isz Ic 

/ check 1 count 

7532 

/5315 


jmp il 

/ repeat inner loop 

7533 

/7325 


7325 

/ 3 

7534 

/1003 


tad um 

/ increment address um 

7535 

/2346 


isz me 

/ check m count 

7536 

/5310 


jmp ol 

/ repeat outer loop 

7537 

/5214 


jmp reg+1 

/ go back to get next set of variables. 

7540 

/0000 




7541 

/0106 

pi 06, 



7542 

/ 02361 

p236, 



7543 

/7300 

awxyz, 


/ address of routine w/xyz 

7544 

/0000 




7545 

/0000 

ac, 


/ a count. lc=ac 

7546 

/0000 

be. 


/ b count. me = be 

7547 

/0000 

cc, 


/ c count 

7550 

/2014 

one, 


/ floating constant — unity 

7551 

/0000 




7552 

/0000 




7553 

/2077 

va', 


/ variable a. 

7554 

/6400 




7555 

/ 0000 




7556 

/2014 

vb, 


/ variable b 

7557 

/0000 




7560 

/0000 



/ variable c 

7561 

/2014 

VC , 


7562 

/0000 




7563 

/0000 




7564 

/2056 

w, 


/ variable w 

7565 

/4000 




7566 

/0000 




7567 

/2035 

X, 


/ variable x 

7570 

/0000 




7571 

/0000 



/ variable y 

7572 

/2014 

y. 


7573 

/0000 




7574 

/0000 



/ variable z 

7575 

/2024 



7576 

70000 




7577 

70000 


um=03 

ul<=05 

xl=06 

ul=07 




QUICK DUMP 


7600 

/6032 qd> , 

kcc 


7601 

/6040 


tfl 


7602 

/7604 


las 


7603 

/ 1267 


tad 

m1 

7604 

/3010 


dca 

10 

7605 

/ 1270 


tad 

m4 

7606 

/3272 


dca 

c4 

7607 

/4255 


jms 

cr 

7610 

/4255 


jms 

cr 

7611 

/4255 


jms 

cr 

7612 

/4255 


jms 

cr 

7613 

/ 1275 


tad 

m4Q 

7614 

/3271 


dca 

c1 

7615 

/ 1265 


tad 

p207 

7616 

/4232 


jms 

tp 

7617 

/ 1010 


tad 

10 

7620 

/7001 


iac 


7621 

/4240 


jms 

pvt 

7622 

/ 1410 


tad 

i 10 

7623 

/4240 


jms 

pw/ 

7624 

/2271 


iS2 

d 

7625 

/5222 


jTip 

• -3 

7626 

/2272 


iS2 

c4 

7627 

/5210 


jmp 

qd+1 0 

7630 

/7402 


hit 


7631 

/5205 


jmp 

qdi+5 

7632 

/7654 

tp. 

0 


7633 

/6046 


tsf 


7634 

/6041 


jmp 

.-1 

7635 

/5234 


tls 


7636 

/3377 


dca 

chr 

7637 

/5632 


jmp 

i tp 

7640 

/7624 

pw. 

0 


7641 

/3346 


dca 

tm 

7642 

/ 1346 


tad 

tm 

7643 

/7002 


bsw 


7644 

/4250 


jms 

tc 

7645 

/ 1346 


tad 

tm 

7646 

/4250 


jms 

tc 

7647 

/5640 


jmp 

i pw 

7650 

/7645 

tc , 

□ 


7651 

/0266 


and 

p77 

7652 

/ 1274 


tad 

p240 

7653 

/4232 


jms 

:tpr 

7654 

/5650 


jmp 

i tc 


/ quick dump 

/ get first address to be dumped 
/ subtract 1 
/ put in index register 

/ c4 + -4 
/ return carriage 

/ nop these two instructions 
/ to single space the output. 

/ set counter at 

/ 32 words per line 

/ start each line with a bell. 

/ and add’ress of first word. 

/ print address of first word 
/ get next word 
/ and type or print word 
/ end of line ? 

/ yes. quit after four lines 

/ one page has been dumped 
/ return to print next page 
/ type character in accumulator 


/ save character 
/ return 
/ print word 
/ save word 

/ type first character 

/ type second character 
/ return 

/ type character 
/ retain low order chara 
/ format for printing 
/ type it 
/ return 
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7655 

/7613 

cr, 

Q 


/ 

carriaga return 

7656 

/7200 


cla 




7657 

/ 1273 


tad 

p21 5 



7660 

/4232 


jms 

tp 



7661 

/ 1264 


tad 

^1212 



7662 

/4232 


jms 

tp 



7663 

/5655 


jmp 

i cr 



7664 

/0212 

p2'12. 





7665 

/0207 

p207. 





7666 

/0077 

P77, 





7667 

/nil 

m1 , 





7670 

/7774 

m4, 





7671 

/7756 

c1 , 





7672 

/7774 

c4. 





7673 

/0215 

P215, 





7674 

/0240 

p24Q, 





7675 

/7740 

m40. 





7676 

/0100 

p1 00, 





7677 

/7700 

m1 00- 





7700 

/6032 

qi, 

kcc 


/ 

quick load 

7701 

/6040 


tfl 




7702 

/ 1353 


tad 

15310 



7703 

/3313 


dca 

gw-2 

/ 

restore last instruction 

7704 

/4315 


jms 

gw 

/ 

get next word from keyboard 

7705 

/5304 


jmp 

.-1 

/ 

return for address if a bell was heard. 

77 06 

/ 1267 


tad 


/ 

subtract 1 from address 

7707 

/3010 


dca 

10 

/ 

put in index register 

77 10 

/4315 


jms 

gw 

/ 

get next word. 

7711 

/5304 


jmp 

.-5 

/ 

return for address if bell was heard. 

7712 

/3410 


dca 

i 10 

/ 

load word in core. 

7713 

/5310 


jmp 

.-3 

/ 

go back for next word 

7714 

/7571 

m207, 





7715 

/771 1 

gw/. 

Q 




7716 

/7240 


cla 

cma 



7717 

/4326 


jms 

rd 

/ 

read first character 

7720 

/7002 


bsw 


/ 

switch to high order position 

7721 

/7421 


mql 


/ 

save it 

7722 

/4326 


jms 

rd 

/ 

get next character 

7723 

/7501 


mqa 


/ 

combine two characters 

7724 

/2315 


isz 

gw 

/ 

since bell was not heard, skip. 

7725 

/5715 


jmp 

i gw 

/ 

return 



7726 

/7723 

rd. 

Q 


/ 

7727 

/4346 


jms 

re 


7730 

/4232 


j?ns 

4* m 


7731 

/1377 


tad 

chr 


7732 

/ 1314 


tad 

m207 

/ 

7733 

/7450 


sna 


/ 

7734 

/5715 


jmp 

i gw 


7735 

/1345 


tad 

m3l 

/ 

7736 

/7510 


spa 



7737 

/5327 


jmp 

rd+1 

/ 

7740 

/ 1277 


tad 

ml OG 


7741 

/7500 


sma 


/ 

7742 

/5327 


jn>p 

rd+1 

/ 

7743 

/ 1276 


tad 

pi 00 

/ 

7744 

/5726 


jmp 

i rd 

/ 

7745 

/7747 

m3l , 




7746 

/7730 

re, 

0 


/ 

7747 

/6031 


ksT 



7750 

/5347 


jmp 

.-1 


7751 

/6036 


krb 



7752 

/5746 


jmp 

i re 


7753 

/5310 

i53lD 

f 



7754 

/0204 





7755 

/0022 




/ 

7756 

75300 




/ 


/ 

/ 

/ 

/ 

/ 


read characters 


is it a bell ? 

yes—(return without a skip) 

no — check is it less than 240. 

yes. ignore and get next char. 

is it greater than 337 
yes, ignore, 
restore character 
return 

get next character 


the first instruction in the 
rim loader will be overwritten 
with a transfer to ql Wfhereupoh 
the instruction will be restored 
and than a halt will be generated 
by bell ^NP: 

bell at end of tape. 
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EXAMPLES 

We will now present two examples to show how the program 


works* 


is a very short problem 


The first, called SMALL TEST , 
involving a logarithmic fit. 


asking for 


Keenan and Keys . John Wiley &, Sons — 1 936, 

In this example, we specified 16 degrees of freedom* 
But, we ran the problem 4 ways to show how the Standard Dev¬ 
iation varies when we specify more degrees of freedom. Results 
are given in the following table; 


called 

BIG 

TEST 

, is a fair sized pr 

oblem 

rees of 

freed 

lorn) . 

In this problem, we 

are 

a with ' 

which 

the 

computer may calculat 

e the 

hen pressure, 

and 

temperature are give 

n .The 

, ^Thermodynamic 

Properties of Steam” 

» lay 


Standard 
A B Deviation 


4 

4 

9 

7 


.000423 
.000238 
.000128 
. 000115 


The computer took less than 40 seconds to 
solve 28 simultaneous equations.It expanded 
the 28 degree polynomial equation in less 
time than it takes for a carriage return. 


From statistics, we learn that , with a normal distribution 
and an infinite number of points, 6B,26% of all cases will be 
within one standard deviation from the mean, and 95,46% will he 
within two standard deviations. In our problem, with 45 points, 
no point exceeded two deviations (.00085) ; but, one point — 

number 18 — was very close. In any case, the accuracy of 
formula is about what you could get by a manual reading of 
Steam Chart. 


the 

the 


SYMBOL 

W 

X 

Z 


We define 

PROPERTY 

ENTHALPY 
AB5. PRES. 
TEMPERATURE 


symbols as follows: 

RANGE 


UNITS 

Thous Btu/lb 
Thous psia. 
Hund. Deg. F 


1,4 - 1,52 
.54- .62 

8.0 - 10.0 


To avoid exponent overflow, it was necessary to scale 
units near unity. For example, in the problem having 28 degs, 
of freedom, the pressure is raised to the twelth power and 
multiplied by temperature to the sixt^gpower. This product 
exceeded the range or our exponent(l0 ) when we tried to 
correlate to a natural scale. 
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REGRESSION ANALYSIS EXAMPLE — SMALL TEST 
Start at address 7342 
DATA 


3A 2D E 
2.5W l.X C 

2.9W 2 .x C 

5.V 4.x C 

I 2 .W 8.x C\ PARAMETERS AND REGRESSION COEFFICIENTS 

7775 7775 7777 7777 7776 0011 

9.130124E-01 -2.918390E-02 3.783139E-01 

Start after halt at address 7265 


RESULTS 


V 

(v - w) 

N 


w 

X 

Y 

Z 


3A 2D E 
2.5W l.X t 

2.491817E+00 
2.500000E+00 
2.9W 2.x C 

2.928662E+00 
2.9 00000E-I-00 
5.V 4 .x C 

4.951067E+00 
5.000000E+00 

12 . W 8.x i 

1.203940E+01 

1.200000E-f01 

3.476840E-02 


-8.182645E-03 
1.000000E-»'00 

2.866203E-02 

2.000000E+00 

-4.893285E-02 

4.000000E+00 

3.940141E-02 

8.000000E+00 


1.000000E+00 

1.000000E+00 

2.000000E+00 

1,000000E*00 

3.000000E+00 

l,000000E+00 

4.000000E+00 

1.000000E+00 


1.000000E+00 

1.000000E+00 

1.000000E-f00 

1.000000E+00 \ 


RMS or Standard Devistion 


Start after halt at address 4420 

DRC = Dump regression Coefficients. QUICK LOAD will put these 
coefficients in your program, 

r ST = Small Test, Put in switch register before start 
at address 4420, See translate table on page 14, 

"••3ST )0’3»ENO?>>B§/t #9 + G 








REGRESSION ANALYSIS EXAMPLE 


BIG TEST 


AA 4B 

i. 4103V 
J.4099W 

j. 4094W 
IV 4 09 W 
.1V4086W 
jV4081W 
JV407TW 
jV4072W 
1V4068V 


DATA 

.54X 8Y 
V55X C 
V56X C 
V57X C 
V58X t 
V59X t 
V60X C 
V61X C 
V62X C 


C 


j.4376W .54X 8.5Y C 

iV4372W V55X f 
JV4368W V56X C 
jV4364W V57X C 
JV4360W V58X C 
1V4356W V59X t 
jV4352W V60X C 
iV4348W V61X t 
IV 4344V V62X L 


i.4646W .54X 9Y C 

iV4643W V55X t 
jV4639W V56X C 
1V4636W V57X C 
jV4632W V58X t 
jV4629W V59X [ 
iV4625W V60X C 
iV4622W V61X t 
1V4618W V62X C 


J.49.15W .54X 9.5Y C 

.1V4912W V55X I 
iV49 09W V56X f 
1V4906W V57X C 
iV4903W V58X t 
iV4900V V59X C 
iV 4897V V60X C 
IV4893V V61X C 
IV 4890V V62X t 


i.5i85V .54X i0Y C 

iV5i82V V55X C 
iV5i79V V56X C 
IV 5176V V57X C 
IV 5173V V58X t 
iV5170V V59X C 
IV5167V V60X C 
IV 5165V V61X I 
1V5162V V62X C 


\ PARAMETERS 

7760 7774 7774 7777 0000 0060. REGRESSION 
2.030484E+00 -7.667827E-01 1.211007E+00 

-2V310170E-01 -2V370390E-01 6V26i075E»01 - 

2V59929 IE-02 5V195604E-02 -9V607082E-02 

-1V128037E-03 -1V461558E-04 - rV507848E-03 


COEFFICIENTS 
9.551607E-01 
9V9 78330E-01 
IVl 19160E-01 
9V50n56E-05 



REGRESSION ANALYSIS EXAMPLE 


BIG TEST 
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RESULTS 



V 

(V - W) 

N 


4A . 4B E 

w 

X 

Y 

Z 


1,.4103W .54X 

■I.4i0078E+00 
. r. 4l0300E-i-00 
1»;4099W _ .55X 
i.409951E-f00 
. lV409900E<f00 
4094W • 56X 

.I.409744E-I-00 
. 1.409400E<i>00 
1,409 W •57X 
I.409 459E<t-00 
. 1V409000E+00 
1» 4086V • 58X 

.I.409096E+00 
. i,'408600E+00 
1,4081V ,59X 

I. 408658E-I-00 

. lV408100E-i-00 
1,4077W • 60X 

X.408144E+00 
.1V407700E+00 
1,4073V ,6IX 

■J.407556E+00 
.1V407200E+00 
1,4068V •62X 

i, 406896E-I-00 
l,406800E+00 

1,4376V •54X 

i,436834E-i'00 
. l,437600E+00 
1,4372V . .S5X 

J, 43662 lE-t-00 

. 1.437200E+00 
1,4368V «56X 

i,436346E-f00 
. r,'436800E-t-00 
1, 4364V. . 57X 
1,43601 1E4>00 
. lV436400E-f 00 
1,'4360V . .58X 
1, 435621E<('00 
. 1V436000E-I'00 
1,4356W ,59X 

■X,435179E+00 
. iV435600E-f00 
1,4352V , 60X 

i,434689E+00 
. 1V435200E-I-00 
1^ 4348V ,6IX 

l,434154E+00 
. 1.434800E+00 
1,4344V •62X 

1.433579E-f00 
r,434400E4>00 


8Y . C 

-2.219826E-04 

5,'400000E>01 

l,000000E+00 

8,'000000E+00 

l,000000E-t>00 

L 

5,067885E-05 

5,500000E-01 

• • - 

S»000000E*00 

8,'000000E+00 

l,000000E+00 

L 

3.439337E-04 

5,600000E-01 

r ■ • 

3,000000E-l-00 

8,000000E^00 

l,000000E-i-00 

*• 

4.587770E-04 

5,700000E-01 

C' 

4.963577E-0.4 

5V800000E-01 

... . 

4,000000E4-00 

8.000000E-i>00 

l,000000E-i'00 

5,000000E-i'00 

8,000000E-i-00 

l*000000E-i-00 

1. . 

5.576760E-04 

5,'900000E-01 

C‘ 

4,439652E-04 

6,000000E-01 

6,000000E-i'00 

8,'000000E-i>00 

l,000000E-f00 

7, 000000E-f00 
8,000000E<<'00 

1,000000E*00 

L 

3,562570E>04 

6,'100000E-01 

8,000000E+00 

8,'000000E-i-00 

l,000000E-(-00 

L 

9.568036E-05 

6’,200000E-01 

9»000000E*00 

8,'000000E'«'00 

'l,000000E->-00 

8,5Y C 
-7,663816E-04 
5,'400000E-01 

l,000000E-f01 

8,500000E-(-00 

l,000000E-f00 

V 

-5,789101E-04 

5,'500000E-01 

f * * 

1, 100000E<)-01 
8,'S00000E-t>00 

l,000000E-l-00 

L 

-4,543066E-04 

B.600000E-01 

f •• 

l,200000E-<-01 

8,'500000E4-00 

1,000000E>00 

-3,890395E-04 

5,*700000E-01 

[■ . 

-3,791600E-04 

5,'800000E-01 

|,*- 

l,300000E-i-01 

8,'500000E-l-00 

l,000000E+00 

1.400000E<»>01 

8V500000E-f00 

l,000000E+00 

t ... 

-4,213601E-04 

5V900000E-01 

r‘* 

l,500000E-i>01 

8','500000E4-00 

l,000000E-t-00 

-5, 114824E-04 

'6V000000E-01 

... 

1,600000E+01 

8','500000E-i-00 

l,000000E4-00 

L 

-6,460995E-04 

6,'100000E-01 

_ . - .... 

l,700000E-t>01 

8,'500000E-i-00 

l,000000E<<-00 

-8,213520E-04 

6,'200000E-01 

1,800000E+01 

8,'500000E-i>00 

l,000000E+00 
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1> 4646V. • 54X 

I.464614£-f00 
.1V464600E+00 
'4643V .SSX 
X.464334E1-00 
.1V464300E+00 
U 4639V. .56X 

I. 464012E+00 

. l.'463900E^00 
1,46 36V .57X 

X.463657E+00 

. l'.'463600E-i'00 
1,4632V .58X 

J. 463274E+00 

. r.463200E-i'00 
1,4629V .59X 

X.462869E-I-00 
. r.‘462900E+00 
1,4625V . .60X 
X.462451E+00 
. 1V462500E+00 
1,4622V *6IX 

X.462024E+00 
. l'<>'462200E-t-00 
1,4618 V •62X 
X.46X597E+00 
l.'46l800E+00 


9Y C 

1.408160E> 05 
5.'400000E>01 
C 

3.378093E>05 

5V500000E-01 

f. 

1.124889E-04 
5.600000E-01 
Z 

5.693734E-05 

5V700000E-01 

c • 

7.365644E-05 

5.'800000E-01 

C' 

•3.077090E-05 

5.'900000E-01 

C" 

•4.942715E-05 

6'.000000E-01 

[■ . 

• 1.758188 E-04 
6 .100000E-01 
C" 

•2.032965E-04 

6.200000E-01 


1.9 00000E-I-01 
9.000000E-i-00 

2.000000E+01 

9.000000E-l'00 

2. 100000E-I-01 
9.000000E+00 

2.200000E-I-01 

9.000000E^00 

2.300000E+01 
9. 000000E-t-00 

2.400000E-f01 

9.'000000E+00 

2.500000E-t-01 

9.000000E-«-00 

2.600000E-I-01 

9.000000E+00 

2. 700000E-«-01 
9.'000000E4>00 


1.000000E+00 

1.000000E<t-00 

1.000000E<«-00 

1.000000E-i-00 

1.000000E4-00 

1.000000E<f00 

l,000000E+00 

1.000000E+00 

1.000000E+00 



1,4915W .54X 

X.492196E+00 
. 1.491 S00E-I-00 
1,'49 12V .55X 

X.49X852E+00 
. r.49l200E-»-00 
1,4909V .56X 

X.491495E-«-00 
. 1.490900E<(-00 
1,49 06W .57X 

X.491134E-<-00 
. r.490600E4>00 
1,49 03V .58X 

X.490779E4-00 
. r.’490300E400 
1,'4900V .59X 

'X.490440E-i'00 
. 1.490000E+00 
1,489 7W •60X 
X.490127E-f00 
. r.*489700E+00 
1, ■4893V .6IX 

X.489850E+00 
. l'.'489300E-i-00 
1,489 0V. .62X 

X.489618E+00 
r.'489000E+00 


9.5Y . Z 

6'.'9 56160E-04 

5V400000E»01 

Z' . 

6.521046E-04 
5'.500000E-01 
C ■ 

5.948693E->04 

5V600000E-01 

C' 

5.337894E-0.4 

5'.700000E‘>01 

Z\ 

4.78 7 594E-04 
5'.800000E>01 

4.396588E-04 
5'.900000E-01 
C ■ 

4.266351E-04 
6.000000E-01 
C" . . 

5.495101E-04 
6.100000E-01 
C ■ . 

6. 181896E-04 
6'.200000E-01 


2.800000E-i>01 

9.'500000E+00 


3.500000E-f01 

9'.500000E-i-00 

3.600000E<l-01 

9'.500000E-l-00 


i,000000E*00 


1.000000E-t'00 

1.000000E-f00 


2.900000E+01 

9'.500000E+00 

3.000000E-i-01 

9'.'500000E4-00 

3.'l00000E-(>01 

9'.S00000E4'00 

3.200000E<i-01 

9.500000E-f00 

3.300000E-t>01 

9'.500000E+00 

3.400000E-I-01 

9'.500000E'^00 


1.000000E‘i-00 

1.000000E-f00 

1.000000E+00 

1.000000E<t-00 

1.000000E+00 

1. 000000E+00 



1,5185V .54X 

X.5i 83541:+00 
. l'.'5l8500E+00 
1.5182V .55X 


10Y Z 

> 1.45629 1E> 04 
5'.400000E-01 
C' 


3.700000E+0X 

r.000000E+0l 


1.000000E+00 
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♦ 


i.5i7940E-i-00 >2.604276E-0A 
. 1V518200E+00 5V500000E-01 

U5I79W .56X C" 
J.517543E+00 -3.566593E-0A 
.1V5I7900E+00 5V600000E-01 

lr5I7$W .57X C • . 
J.Si7l79E+00 -4,810770E-0^ 
. lV5t7600E<i-00 '5V700000E-01 
i,5i.73w .sex C V 
X.S16860E+00 -4.399270E-0.4 
. IV5l7300E-i^00 5V800000E-01 

U51.70W .59X C* 

I.SJ6600E+00 -3.997237E-S4 
. 1 VS 17000E<i>00 5V9 00000E- 01 

1V5IS7W. .60X C" 
I.S16413E-I>00 -2.867729E-04 
. rVSl6700E>f00 6V000000E-01 

1,5165V. .SIX C . 
I.Sl6312E<t-00 >U876354E>04 
. rV5l6500E4'00 6V100000E>01 

1V5162W.. .62X r;. 

'i.5i63UE-i>00 UU4309E-04 

1V516200EV00 6V200000E-01 

^ RMS or DRC 

Standard Deviation 

4.226363E-04 


3,800000E<t-0i 

rV000000Et0l 

!• 000000E’i>00 

a»900000E*01 

1V000000E+01 

l«000000E-i-00 

4«000000E+0i 

1V000000E*I-01 

1.000000E-I-00 

4*V00000E+0i 

1V000000E+01 

l,000000E*00 

4 b 200000E+ 01 
1V000000E401 

1*000000E*00 

4.300000E-I-01 

1V000000E+01 

1.000000E+00 

4.400000E-i-0j 

1V000000E+01 

1.000000E+00 

4.500000E+0i 

rV000000Et01 

1.000000E-f00 


= Dump Regression Coefficients. 
QUICK LOAD will put these 
coefficients in your program^. 


— BT = Big Test. Put in switch register before 
start at 4420. See translate table on 
page 14. 

P04#YN|P&<ECV0>V!)00*I"M+QW9!VGOWE»+X0X D2?P*tX->/>IW 
BC/FIQC90N>«8\/0*< •KN\GM46NDR9M(Nt+l>«-.>.®>+ 








% 


50 


ADDRESS SUl-lMARY 




27 

BIT FLOATING POINT 

' PACKAGE 






5200-6777 





PROG. 

FIX 

ANTE 

RINT 

FIN 

LOG 

FOUT 

EDIT 


ADDR. 

5234 

5271 

6502 

6600 

6661 

6200 

5200 








--HALTS 



PROG. 

TYPE 

C. R. 

POA 

FPNT 

EXP. 

EXP. 

DIV. 


ADDR. 

6305 

5066 

5110 

5400 

OVER 

L'N DER 

BY 

NEG. 

INST. 

4475 

4474 

4473 

4476 

FLOW 

FLOW 

ZERO 

LOG 






5525 

6046 

6064 

7402 




TRACE : 

ROUTINE 








5034-' 

-5177 







-TRACE- 







LINK 

LIVE 

DEAD 






INSERT 

5134 

5477 

2200 

LOW 

HIGH 




e ADDR. 

0077 

5404 

5404 

5050 

5051 







TEf'IPORARY STORAGE ON 

PAGE ZERO 

PROG. 


z 



LOG 

VAR. 

Y 

F INT CNT 

CH 

ANTE 

ADDR. 

35-37 

40-42 43-45 46 

47 



FPNT 

--SIGN- 

— 

-EXPONENT- 



VAR. 

OP 

AC 

OP 

AC 

AC 3- 1 

OP 3-1 

ADDR. 

50 

51 

52 

61 

53-55 

56-60 





